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ABSTRACT 

The shutdown of star formation in galaxies is generally termed 'quenching'. Quenching may occur through 
a variety of processes, e.g., AGN feedback, stellar feedback, or the shock heating of gas in the dark matter halo. 
However, which mechanism(s) is, in fact, responsible for quenching is still in question. This paper addresses 
quenching by searching for traces of possible quenching processes through their effects on galaxy structural 
parameters such as stellar mass (M,), M,/re, surface stellar mass density (~ M^/r^), and Sersic index (n). We 
analyze the rest-frame U-B color correlations versus these structural parameters using a sample of galaxies in 
the redshift range 0.5 < z < 0.8 from the DEEP2/AEGIS survey. In addition to global radii, stellar masses, and 
Sersic parameters, we also use 'bulge' and 'disk' photometric measurements from GIM2D fits to HST/ ACS V 
and / images. We assess the tightness of the color relationships by measuring their 'overlap regions', defined 
as the area in color-parameter space in which red and blue galaxies overlap; the parameter that minimizes these 
overlap regions is considered to be the most effective color discriminator We find that Sersic index (n) has 
the smallest overlap region among all tested parameters and resembles a step-function with a threshold value 
of n = 2.3. There exists, however, a significant population of outliers with blue colors yet high n values that 
seem to contradict this behavior; they make up « 40% of n > 2.3 galaxies. We hypothesize that their Sersic 
values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central 
surface stellar mass density, I^*]^^^, as an alternative to Sersic index. Not only does 1^*^^^. correct the outliers, it 
also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically 
linked with quenching. Furthermore, at z ~ 0.65, the majority of the blue cloud galaxies cannot simply fade 
onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses 
of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers 
of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates 
are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in 
later, associated with the central stellar mass build-up. The mass build-up is naturally explained by any non- 
axisymmetric features in the potential, such as those induced by mergers and/or disk instabilities. However, the 
identity of the second quenching agent is still unknown. We have placed our data catalog on line. 

Subject headings: galaxies: bulges — galaxies: formation — galaxies: evolution — galaxies: structure — 
galaxies: fundamental parameters 



|ec2250@gmail.com] 

' Department of Astronomy and Astrophysics, 1156 High Street, Uni- 
versity of California, Santa Cruz, CA 95064 

^ UCO/Lick Observatory; University of California, 1156 High Street, 
Santa Cruz, CA 95064 

Department of Physics and Astronomy, Elliot Building, 3800 Finnerty 
Road, University of Victoria, Victoria BC, V8P 5C2, Canada 

'' Herzberg Institute of Astrophysics, National Research Council of 
Canada, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden State, Cam- 
bridge, MA 02138, USA 

* Steward Observatory, University of Arizona, 933 North Cherry Av- 
enue, Tucson, AZ 85721, USA 

' National Optical Astronomical Observatories, 950 N. Cherry Avenue, 
Tucson, AZ 85719, USA 

^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, 
MD 21218 

' Department of Physics and Astronomy, University of Pittsburgh, 3941 
O'Hara Street, Pittsburgh, PA 15260, USA 

'" Department of Physics, Center for Astrophysics and Space Sciences, 
University of California at San Diego, 9500 Oilman Dr. La Jolla, San 
Diego, CA 92093 

' ' Kavli Institute for the Physics and Mathematics of the Universe, Uni- 
versity of Tokyo, Kashiwa, 277-8582, Japan 



'- University of Nottingham, School of Physics and Astronomy, Not- 
tingham NG7 2RD, UK 

''' Department of Physics and Astronomy, University of California, 
Riverside. CA 92521, USA 

'* Center for Galaxy Evolution, Department of Physics and Astronomy, 
University of California, Irvine, 4129 Frederick Reines Hall, Irvine, CA 
92697. USA 

'''Max Planck Institute for extraterrestrial Physics, PO Box 1312, 
Giessenbachstrasse, 85741 Garching, Germany 

'* Racah Institute of Physics, The Hebrew University, Jerusalem 91904, 
Israel 

'^ Astrophysics Science Division, Goddard Space Flight Center, Code 
665, Greenbelt, MD 20771 

"* Department of Astronomy, University of Michigan, 500 Church St., 
Ann Ai-bor, Ml 48109, USA 

" Steward Observatory, 933 N. Cherry St., University of Arizona, Tuc- 
son, AZ 85721, USA 

Institute of Astronomy and Astrophysics, Academia Sinica, Taipei 
106, Taiwan 

Department of Astronomy, University of California, Berkeley, CA 
94720 

NASA Postdoctoral Program Fellow 
t Hubble Fellow 



2 



Cheung et al. 



1. INTRODUCTION 



With the advent of large galaxy surveys, the color 
bimodality of the galaxy population has become well- 
characterized jLin et all 1 19991: IStrateva et alJl2001 |: llm et al " 



200l 



Blanton et all 120031: iKau ffmann et al. 2003r lBell et'al 



2004bl) . Galaxy counts back in time revealed that the num- 
ber of red galaxies has at least doubled since z ^ I while 
the number of blue galaxies has remained relatively con- 
stant (Bell et al. 2004b; Bundv et al. 2006; Faber et al 2007; 
Arnoutsetal. 2007; Brown et al. 2007; llbert et all boia 
Dommguez et al..,201 l:.Gon9alves et al.,,2012.) . A natural in- 



terpretation is that galaxies evolve from blue to red with time, 
i.e., from star-forming to 'quenched'. Later measurements of 
star formation rates confirmed that blue galaxies create stars 
at a high rate wh ile red galaxies show little t o no star forma- 
tion (ISahm etaii,2005 , 2007; Noeske et al...2()07HZheng et alJ 



l2007h . Moreover, star formation rates in blue galaxies cor- 
relate well with stellar mass, forming the 'Main Sequence' 
of star formation. In non-dusty red galaxies, however, star 
formation is gene rally much lower than that of blue ones 
(ISalim et alJl2005l) . This abrupt jump in star formation rate 
across colors motivates the search for a quenching process. 
For simplicity, we define quenching to be a process that per- 
manently turns a blue star-forming galaxy into a red non-star- 
forming one PI 

Many quenching mechanisms have been proposed, but 
they can generally be categorized into two classes. The 
first class is internal processes; these act to either expel 
the gas already in a galaxy or render it inert to star for- 
mation. Examples of internal processes include feedback 
from starbursts and active galactic nuclei (AGN), both of 
which may be triggered by mergers. They act to heat the 
surroundin g gas and/or drive winds out of the galaxy (e.g.. 



Sanders et a l. 1988; Springel et al. 2005; Murrav et al. 2005|; 
Cox et all 12008^ .Ciotti et al., ,2009; Alexander et al., ,20 id) . 



Another ex ample of an inter nal process is morphological 
quenching ("Martig e t al.ll2009h . In this model, the presence 
of a dominant bulge stabilizes the gaseous disk against gravi- 
tational instabilities needed for star formation. 

The second class contains external processes, which we 
define as acting to prevent gas from accreting onto a 
galaxy in the first pla ce. The main external process is 
halo mass quenching dSiUd '1977; Rees & Ostriker 1977^ 
Blumenthal et al. 1984; Birnboim & Dekel 2003; Kere s et all 
2005; Dekel & Birnboim 2006; Cattaneo et al. 2006); this 
posits that dark matter halos above a critical halo mass es- 
tablish virial shocks that stop the flow of cold gas onto 
their central gal axies. Additional examples are AGN 'radio 
mode' feedback (Croton et al. 2006) and gravitational heating 
liChochfar & Ostriker 2008; Birnboim & Dekel 2011), both 
of which can be considered as variants of halo mass quench- 
ing since both mechanisms require massive halos. 

According to our definition, mergers do not qualify as an 
external process since they act to exhaust and/or remove ex- 
isting gas already wit hin a galaxy. By the same token, ram 
pressure stripping (Gun n & GottI 11972) is also not consid- 
ered an external process since it strips gas from a galaxy. 
Furthermore, this paper only concentrates on quenching pro- 
cesses that affect the central galaxy of a halo. According to 

Rejuvenation of star formation in quiescent spheroids through 
gas and/or satellite infall has been proposed t o explain the observed 
blue sp heroids seen in various works (e.g., IKannappan et al.l 120091 : 
ISchawins ki et al. 20W). In this paper, we do not consider this process. 



iGerke et akl (12005 1). who used a sample of DEEP2 galaxies 
(~ 25% of the total DEEP2 sample), only - 32% of DEEP2 
galaxies are in groups, meaning ^68% of these galaxies are 
in the field. These field galaxies would be centrals and ad- 
ditionally, since each group contains one central, the per- 
centage of centrals in this DEEP2 sample is at least ^68%. 
Assuming this sample is representative of the entire DEEP2 
dataset, we can conclude that most of our galaxies are cen- 
trals. Thus we will not consider mechanis ms that affect satel- 
lites, i.e., st rangulation and harassment (iLarson et al.lll980l 
iMooreet al . 1996). 

These quenching processes may imprint themselves on the 
structure of a galaxy, e.g., major majors can create highly 
concentrated galaxies. The prospect of detecting quenching 
mechanisms at work via observable changes in structural pa- 
rameters has motivated many previous works. One of the first 
parame ters explored was luminosity. Using an early SDSS 
sample. [Strateva et al . (2001) found that galaxies are bimodal 
in color, i.e., galaxies generally lie within the red sequence or 
the blue cloud. However, while galaxies are well separated in 
color, they overlap over almost the entire range of luminos- 
ity, indicating that luminosity is not the main driver of galaxy 
color 

Later, stellar mass was explored; hereafter, mass refers to 
stellar mass un less otherwise stated. F or a sample of local 
SDSS galaxies. IKauffmann et al.l ( 12003b found that the corre- 
lations between the star formation history indicators Dn(4000) 
and HSa (which can also be thought of as a proxy for galaxy 
color) and mass are significantly better than that of the g-band 
luminosity. They further found that galaxies divide into two 
distinct families at a stellar mass of 3 x lO'*' M©. 

Recently, additional structural parameters have b e en in- 
troduced. Using an SDSS sample, Kauffma nn et al.l (120061) 
found that the galaxy surface mass density (^ M^/r^) pro- 
duced an even sharper division in s pecific star formation rate 
(SSFR) than stell ar mass (see also iBrinchmann et al.l 120041 ; 
iMaier etal] 12009 ). They suggested that high surface stellar 
mass density is connected to the creation of a bulge and the 
quenching of a galaxy. 

Franx et al. (2008) intercompared several of the aforemen- 
tioned color-parameter correlations in the redshift ran ge < 
z < 3 using data from the FIREWORK S catalo g (.Wuvts et alj 
120081) . Confirming IKauffmann elal] (I2006l) 's result, fliey 
showed that surface mass density better separates red and blue 
galaxies than stellar mass alone. Franx et al. (2008) also ex- 
amined a second structural parameter, the "inferred velocity 
dispersion" (~ M*/re), and found that the inferred velocity 
dispersion also better distinguishes red and blue galaxies than 
mass. 

Besides these struct ural parameters, Ser sic in dex (n) has 
also been explored. iDriver et'al] (120061) and lAllen et akl 
(2006) observed a clear bimodal distribution in both the rest- 
frame u- r color and n in the Millennium Gala xy Catalog. A 
similar tren d with SDSS galaxies was seen by 'B lanton et alj 
(I2OOI and iSchiminovich et all dlOOl . Bell (200"1) showed 
that n is an even better color discriminator than surface mass 
density. However, several outliers were noted, and he con- 
cluded that high n is a nec essary (but not suffic ient) condition 
for qu iescence. Recentlv. I Wuvts et al.l (120111) and iBell et al.l 
(I2OI2I) found that the correlation between quenching and n 
was in place s ince at least z ^ 2.5. 

A study by iMendez et alj ( 1201 II) supports the implications 
of the relationship between Sersic index and quiescence. Us- 
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ing a sample of DEEP2/ AEGIS galaxies at 0.4 < z < 1 .2, they 
compared the morphological parameters (CAS, G/M20, and 
B/T) of galaxies in the green valley - galaxies with colors 
that lie between the blue and red peak in the color bimodality 
- to those in the blue cloud and red sequence. They found that 
most green valley galaxies are still disks but are building up 
their central bulge, in that they have higher concentrations and 
higher B/r ratios than blue galaxies and less than red galax- 
ies. In other words, they found that the bulges of galaxies are 
being created or augmented in the evolution of a galaxy from 
the blue cloud, through the green valley, and finally onto the 
red sequence. 

A recent study bv IWake etall (l2012ah adds SDSS central 
velocity dispersion to the list of previously considered struc- 
tural parameters. It is also the first study to compare the ef- 
ficacy of Sersic indices head to head versus other variables. 
They find that central velocity dispersion leaves the weak- 
est residual color trends with other parameters and conclude 
quenching correlates most strongly with central velocity dis- 
persion. 

The dependence of quiescence on halo properties has thus 
far been measured only statistically, by looking at the prob- 
ability that a galaxy is que nched as a functio n of some mass 
and/or surrounding density. iPeng et all ( 120 IQ) found that just 
two processes - "stellar mass quenching", which correlates 
directly with galaxy stellar mass, and "environmental quench- 
ing", which correlates directly with local environmental den- 
sity - can accurately describe the quench ing probabilities of 
SDSS galaxies. A later paper (Peng et al. 201 1) divided cen- 
trals from satellites and found that central quenching - of rele- 
vance here - had no environment de pendence but relate d only 
to stellar mass. A simila r study by IWoo et alj (120121) intro- 
duced halo mass, which iPeng et al.l (120101) had not consid- 
ered, and found that central quenching correlated better with 
halo mass than with stellar mass. However, it is important to 
note that, regardless of whether halo mass is better than stellar 
mass, it is clearly not as predictive as structural variables such 
as Sersic index or central velocity dispersion. We expound 
on this statement in the discussion of this paper, but a cur- 
sory examination of the SSFR as a function of halo mass from 
iConrov & Wech sler (2009) (Fig. 8) shows that star formation 
only gradually changes as a function of halo mass. Whereas 
the plots of color as a function of Sersic index and central ve- 
locity dispersion from I Wake et al.l (12012 a) (Fig. 1) show that 
color changes quite sharply as a function of both these param- 
eters. Thus a central challenge has emerged for the halo mass 
quenching picture, namely, why do galaxy structural param- 
eters predict the outcome of halo mass quenching better than 
halo mass itself does? We return to this question below. 

While correlations do not necessarily imply causality, they 
are strong hints, all of which has led to a rather complicated 
picture of galaxy evolution. Quenching may well involve a 
mix of complex processes that are likely to be dependent on 
several parameters that are themselves correlated. However, 
several themes emerge from the results discussed. The con- 
ditions of the bulge and perhap s the very center o f the galaxy 
appear to be important. Indeed. iKauffmann et al.l (2006) sug- 
gested that bulge-building is the underlying cause of their cor- 
relation between color and surface mass density. And several 
authors, cited above, concluded that high n is necessary to 
quench a galaxy, providing further evidence for bulges. More- 
over, since a hallmark of bulges is high central density, it 
is notable that ' Wake et aTl ( l2012ah find that central velocity 
dispersion is the single most correlated parameter of all with 



galaxy color And finally, since bulges and hi gh central densi- 
ties a r e closely associated with black holes dMagorrian et alj 
119981: iGebhardt et alJlIOOOl) . it is tempting to conclude that 
this mounting chain of evidence is simply a "smoking gun" 
pointing to AGN feedback. 

In total, these works suggest that internal processes, and 
specifically central processes, are responsible for shutting 
down star formation. As noted, this poses a problem at 
first sight for halo mass quenching, since halo properties are 
seen to correlate more weakly with quenching than do vari- 
ables such as Sersic index and central velocity dispersion. 
However, a key element in the ha lo picture is radio mode , 
which depends on AGN feedback (Dekel & Birnboiriil 120061 : 
ICroton et al.ll2006h . and thus possibly on internal/central con- 
ditions. Perhaps it will be possible to link these various pro- 
cesses in a plausible causative chain that explains all of the 
data. We return to this possibility in the discussion section. 

In this paper, we build on previous works and consider the 
possibility of multiple physical processes acting together in 
concert to quench star formation in galaxies. Whereas most 
works have only explored global structural parameters, we ex- 
plore both global and central structural parameters. Our data 
set is AEGIS galaxies possessing HST/ACS imaging, simi- 
lar to fMendez et all (1201 1ft but over a narrower redshift range, 
0.5 < z < 0.8. Like them, we use color as a proxy for quench- 
ing and structural parameters derived from the same GIM2D 
fits. However, we focus on different structural parameters and, 
importantly, convert luminous quantities of subcomponents to 
stellar mass using color-derived M/L ratios. 

Our ultimate goal is to identify that parameter, or combi- 
nation of parameters, that seems to be the best discriminant 
between star-forming and quenched galaxies. Having found 
that combination, we compare its effica cy (or sharpne ss) to 
studies focusing on halo parameters (e.g. JWoo et alJl2()12l) in 
order to assess whether the primary driver of quenching is 
conditions that exist inside a galaxy or outside it. 

A major result of this paper is that the Sersic index, «, dis- 
plays the sharpest break between star-forming and quenched 
galaxies, i.e., it looks most like a quenching threshold. How- 
ever, n does not really distinguish red and blue galaxies all 
that well - ss 40% of AEGIS galaxies with high n have blue 
colors. Suspecting contamination from starbursts, AGN, or 
errors of measurement, we introduce a novel parameter that 
is closely related to n but is more robustly measured, namely, 
central surface mass density, Under this parameter, we 

find that the number of outliers is dramatically reduced, im- 
plying that the innermost structure of galaxies may be most 
fundamentally related to quenching. Moreover, using stellar 
mass measurements of the bulge and inner 1 kpc region of 
galaxies, we show that at z ^ 0.65, most blue galaxies cannot 
simply fade onto the red sequence; they must instead undergo 
a significant restructuring of their innermost stellar density 
profiles en route to quenching. 

These results are compared to various theoretical models. 
The first major conclusion is that the quenching sharpness 
found with our new parameter, central surface mass density, 
far exceeds that found with halo mass, highlighting a major 
tension with the halo quenching picture. Looking at alterna- 
tive theories, we find several striking points of agreement with 
the major-merger picture, but also some important caveats. 
These concerns suggest that bulge-building may just proceed 
quite naturally because galaxies at these redshifts are not yet 
very axisymmetric and non-central torques are constantly be- 
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ing generated. Finally, the very close connection between 
quenching state and central conditions that we find in this pa- 
per looks like a "smoking gun" for AGN feedback, yet not all 
aspects of the data are fuUyexplained by that model either. 

Finally, we place onlinq^ one of the most comprehensive 
datasets available, comprised of 11, 223 galaxies at 0.2 < z < 
1.2, with a mean redshift of z ^ .75. One powerful aspect of 
this dataset is the use of multi-color //5r/ ACS V- and /-band 
imaging, which allows the accurate conversion of light to stel- 
lar mass. It also includes GIM2D bulge-disk decompositions 
(|Simard et al. 2002), which provide photometric and struc- 
tural measurements of the bulges and disks separately; these 
intermediate redshift galaxy decompositions are only possi- 
ble thanks to the high resolution HST imaging. Additionally, 
stellar masses are derived for the subcomponents using their 
y,/ colors. 

This paper is organized as follows: §2 describes our data 
and derivation of the analyzed quantities. In §3, we explain 
our sample selection criteria and discuss sample complete- 
ness. §4 presents our main results - the correlations between 
structural parameters and color In §5, we compare our results 
with several theoretical models and present our two-stage sce- 
nario of galaxy evolution. Finally, we list our conclusions in 
§6. A cosmology with //q = 70 km s"' Mpc"', Vl„, = 0.30 and 
JIa = 0.70 is used throughout this paper All magnitudes are 
on the AB system. 

2. DATA 

We start with a description of all the main sources of data 
used in this paper, which come from AEGIS, and then discuss 
the sample sele ction in ^ For a n overview of the AEGIS 
data, please see lDavis et al.l (I2007h . 



2.1. CFHT BRI Photometric Catalog 

The fist photometric data are we use from the Canada- 
France-Hawaii Telescope (CFHT) BRI imaging catalog. The 
CFHT 12k camera has a 12,288 x 8192 pixel CCD mosaic 
array and a plate scale of 0.207" per pixel, providing a field 
of view of 0.70° x 0.47°. Five separate fields, with one to five 
distinct CFHT 12k pointings per field, were observed from 
1999 September to 2000 October The integration time for 
each points was ^ 1 hour in B and R and for ^ 2 hours in /, 
broken down to individual exposures of 600 s. The data are 
complete to ^ 2 5.25 in B, 24.75 in R, and ^ 24.25 in / (see 
ICoil et al.ll2004l for more details). 

These BRI magnitude s were used with ^-correct v4.2 
jBlanton & Roweisll2007l) to obtain the rest-frame color {U - 
B) and absolute magnitudes {Mb) used throughout this paper. 

2.2. HST ACS V+I Imaging and SExtractor Photometry 

The main photometric catalog from which the sample was 
selected is based on HST/ ACS images taken as part of the 
AEGIS survey (Davis et al. 2007) under program GO-10134 
(PI: M. Davis). The exposures were taken between 2004 June 
and 2005 March over 63 tiles covering an area approximately 
10.1'x 70.5' in size. Each tile was observed for a single or- 
bit in F606W (V) and F814W (/) using a four-point dither 
pattern. These pointings were combined with the STSDAS 
Multidrizzle package using a square kernel. The final images 
have a pixel scale of 0.03" per pixel and a point-spread func- 
tion (PSF) of 0.12" FWHM. The 5a hmiting magnitudes for 

^^ |http : / /people . ucsc . edu/~echeungl/data . html| 



a point source are V = 28.14 and / = 27.52 within a circular 
aperture of radius 0. 12"(^ 50-pixel area). For an extended ob- 
ject, the 5(7 limiting magnitudes are V = 26.23 and / = 25.61 
for a circular aperture of radius 0.3"(~ 314 pixel area). 

SExtractor ( B ertin & Arnout s 1996) is used to detect ob- 
jects in summed ACS V+I images and to construct initial 
galaxy segmentation maps. A detection threshold of 1.5cr and 
50 pixe ls is chosen. These detection maps and the ACS zero 
points dSirianni et alj|2005l) were applied to each band sepa- 
rately to create the ACS photometric catalogs. We selected all 
nonstellar objects with SExtractor CLASS_STAR < 0.9 and 
/ < 25.0 that did not lie within 50 pixels of a tile edge for our 
automated morphology analysis, covering an effective area of 
710.9 arcmin^ in the ACS images (see lLotzetani2008al for 
more details). 

This high resolution catalog was used to generate the galaxy 
sample comprising the GIM2D bulge-ndisk catalogs. 

2.3. GIM2D 

Structural parameters of the HST/ACS imaged galaxies 
were measured using GIM2D, a 2D bulge-ndisk decomposi- 
tion program (Simardet al. 2002). Three separate fits were 
made: a single Sersic fit with floating «, a bulge + disk fit with 
ribuige = 2 and niHsk = 1, and a bulge + disk fit with ntmigf, = 4 and 
ndisk = 1 The three fits were done simultaneously using both 
the V and I HST/A CS images according to the procedure in 
ISimard et alj (120021) . The bulge surface brightness profile is 
parameterized by: 



I](r) = Eeexp{^[(r/re)i/"-l]} 



(1) 



as given by ( ISersicll968h . Here, the parameter k is set equal to 
1.9992n-0.3271, so remains the projected radius enclosing 
half the light (C apaccioli 1989) . The disk profile is a simple 
exponential: 



S(r) = Eoexp(-r/rd), 



(2) 



where Eq is the face-on central surface brightness and is the 
semimajor axis scale length. For the single Sersic fit, Eqn. [T] 
is used to fit the whole galaxy. 

GIM2D also measures concentration, which, unlike the 
SDSS definition, is defined as the ratio of the inner and 
out er isophote fluxes of n ormalized radii a and 1; we fol- 
low l^r^ametaD (119941) and use a = 0.3. Additionally, the 
GIM2D models produce galaxy, bulge, and disk V,I magni- 
tudes - these are the primary magnitudes used throughout this 
paper. 

Throughout this paper, most of our results utilize the sin- 
gle Sersic fits (TableO. When examining the bulge and disk 
properties, we use the best-fitting, two-component decompo- 
sition, i.e., either «/,u/g(, = 4 or nia/ge = 2 (Tables |4] and |5]), for 
each galaxy as indicated by x^- We only use bulge measure- 
ments of Z?/r > 0.1 galaxies, as the low signal-to-noise makes 
measurements of systems with B/T < 0. 1 uncertain. Our sub- 
component sample with GIM2D measurements of the bulge 
and disk separately is comprised of « 60% from the ntuige = 2 
and « 40% from the «/,„/g(, = 4 fits. Comparing the nhuige = 4 fit 
to the nhuige = 2 fit shows a median offset of log bulge to be 
« 0.10 dex with a dispersion of « 0.23 dex while log re.buige 
has a median offset of « 0. 15 dex with a dispersion of « 0.28 
dex; both parameters are offset toward higher values in the 
nbulge = 4 fit. 
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2.4. DEEP2 + DEEPS Redshift Survey 

Spectroscopic redshifts were measured in the DEEP2 
redshift survey using the DEIMOS spectr ograph d Faber et alj 
2003) on the Ke ck II telescope dDavis et al. 2003r 
Newman et alj 120121) . Targets were selected for DEEP2 
s pect roscopy from the CFHT BRI imaging described in 
32.11 Most of DEEP2 used the BRI photometry to screen 
out low-redshift galaxies, but this screening was not ap- 
plied in the AEGIS region, and so the resulting sample is 
representative from z = to z ~ 1 .4. Eligible targets must 
have 18.5 ^ R < 24.1 and surface brightness br ighter than 
IJ,R < 26.5 (iDavis et al.ll200llNewman et al.ll20T2l) . 

Additional spectroscopic redshifts are available in the re- 
cently completed DEEP3 redshift survey (Cooper et al. 20111 
120^121) . This survey shares many of the same characteristics 
of DEEP2, i.e., they both use the DEIMOS spectrograph and 
were both preselected using CFHT BRI photometry. How- 
ever, while DEEP2 used a 1200 line/mm grating in DEIMOS, 
DEEP3 employed a 600 line/mm grating, resulting in spectra 
of lower resolution. The quality of the redshifts, however, are 
unaffected. 

Taking only spectroscopic redshifts with quality code of 
2 = 3 or 2 = 4 and cross matching it to the HST/ACS cat- 
alog yields a sample of 6310 galaxies; these galaxies make up 
the spectroscopic sample. 

2.5. Photometric Redshifts 

The DEEP2H-DEEP3 survey is approximately 65% com- 
plete to R = 24.l in AEGIS (iNewman et all l2012h . For 
those galaxies without spectroscopic z and to extend the sam- 
ple to fainter limits, we utilized photometric redshifts (J. 
Huang et al., in prep) de rived from the Artificial N eural Net- 
works method (ANNz, I CoUister & Lahavl 120041) using the 
multi-wavelength AEGIS photometry that includes 12 unique 
bands in the wavelength range from m to 8 /im, with deep 
Spitze rflRAC photometry (Davis et al. 2007; Barmby et aQ 
120081; Zheng et al., in prep) as the base. This sample was 
3.6 fim selected (/s.e > 2/iJy) with a color cut to isolate 
z < 1.5 galaxies. The redshift catalog is complete down to 
log M^,/Mq = 9.5 for 0.4 < z < 1.2, and the rms accuracy 
is Az/(l+z) = 0.025. Cross matching this sample to the 
HST/ACS catalog that does not have a quality spectroscopic 
redshift yields 4913 galaxies; these galaxies make up the pho- 
tometric sample. The total number of galaxies in our initial 
sample, consisting of the spectroscopic and photometric sam- 
ple, is 1 1,223 

2.6. Rest-frame Absolute B Magnitudes and U — B Colors 

Rest-frame absolute Mb magnitudes and U -B colors are 
needed for both integrated galaxies and for bulge and disk 
subcomponents separately. For galaxies, these quantities are 
obtained t hrough ^ -correct v4.2 (A:-corrected down to z = 0; 
iBlanton & RoweisI 12007) with CFHT BRI photometry and 
redshift as inputs. 

For bulges and disks, however, CFHT BRI photometry is 
not available, but there is HST/ACS V and / photometry mod- 
eled by GIM2D. In order to be consistent with the galaxy val- 
ues, we derive a calibration for Mb and U -B from V, I, and 
redshift. We use the galaxy rest-frame magnitudes from k- 
correct as fiducial values to derive this calibration, which was 
then used to calculate Mb and U -B for the subcomponents 
sample. 

The functional form we use for Mb is (iGebhardt et alj 



I2OO3IPI : 

MB = hi4-DM(n,„,nA,nK)+KjB, (3) 

where DM is the distance modulus for the adopted cosmology 
and 

/T/B = 1 .490 - 1 8 .266z + 94.056z^ 

- 229 .782z^ + 294.74 Iz"* - 1 89 . 892z^ + 48 .034z'^ 

+ (2.233 - 5 .448z+ 3 . 1 87z^- 0.082z^)(V - /) 

+ (0.592 - 0.540Z- 0.036z^)(y - if. (4) 

The functional form for U-B is (IGebhardt et al.ll2003h : 

f/-B = -0.882+ 16.627z-84.798z^ 

+ 212.831z^-286.211z'*+195.256z^-52.584z^ 
+ (0.492 + 0.380z+0.415z^-0.493z^)(y-/) 
+(0.751 -1.609z+0.739z^)(y-/)^ (5) 

Fig. [T^ and[T]5 compares the values of Mb and U-B from 
Eqn. [3]|5] to those derived from fc-correct. The relations are 
nicely Unear with (t(Mb) = 0.215 mag and aiU-B) = 0.096 
mag. We use these relations to compute Mb and U-B for 
the subcomponent sample. We also use these equations for 
the galaxies in our sample that have ill-measured CFHT BRI 
measurements, characterized by large errors (« 7%). 

2.7. Stellar Masses 

Stellar masses for most of our sample are available from 
J. Huang et al., (in prep.). Using a Salpeter initial mass 
function (IMF), the multi-wavelength AEGIS photometry 
(with deep Spitzer/IRAC photometry as the base (! Davis et"al] 
12007; Barmby et al. 2008; Zheng et al., in prep; J. Huang 
et al., in prep) were fit to a grid of synthetic SEDs from 
,Bruzual & Chariot! (120031 ). assuming solar metallicity. These 
synthetic SEDs span a range of ages, dust content, and ex- 
ponentially declining star formation histories. The typical 
widths of the stellar mass probability distributions are 0.1- 
0.2 dex. 

To obtain stellar masses for the subcomponent samples we 
utilize the well-known correlati on between mass-to-lig ht ratio 
(M,/L) and optical colors (e.g.. lBell & de Jong|l2001h . To ac- 
count for our large redshift range, we add a redshift- dependent 
term t o the relationship, similar to the approach in iLin et al.l 
( 120071) and[Weiner et al. (2009). Options are to use either rest- 
frame U-B (from ^-correct) or observed V -I. To aid our 
choice, we make fits using both colors and compare them to 
the M^/Lb values in Fig. |2] The left panels display M^JLb 
vs. observed Y -I while the right panels plot M^,/Lb vs. rest- 
frame U-B. Each row represents a different redshift range. 
Overplotted in each panel is a red dashed curve that repre- 
sents our fit. The fit for M^,/Lb as a function of V-/ is a 
better match than the fit to U-B, especially at high redshifts. 
The final adopted expression for M^/Lb is: 

log M*/Lb = -0.340-2.593z+1.195z^ 

1.908(y-/)-0.432(y-/)2 (6) 

We use lGebhardt et al] <2003l) 's form but fit for our own coefficients 
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Figure 1. We plot the residuals of rest-frame absolute B-band magnitude 
Mb, rest-frame U —B, and stellar mass M, against rest-frame U —B from k- 
correct. Panel a shows the difference between Mb computed from /;-correct 
v4.2 (Blanton & Roweis 2007) using CFHT BRI and redshift and Mb com- 
puted from Eqn.[3]and|4] Panel b plots the residuals of U-B from /:-coiTect 
to U — B from Eqn.|5] Panel c shows the residuals of M, from J. Huang et 
al., (in prep) and M» obtained using Eqn.|3]and \6\ From the one-to-one red 
dashed line in each panel, it is clear that values derived from our fits are con- 
sistent with the fiducial values. The dispersion cr is displayed on the top of 
each panel. 

Together with the absolute B magnitudes measured from 
Eqns. |3] and |4] we are now in a position to calculate stel- 
lar masses for any object with a measured V,I, and redshift. 
These calibrated fits are used to obtain values for the sub- 
component samples. For the galaxies in our sample that do not 
have stellar masses from J. Huang et al., (in prep) (w 10%), 
their M, are also obtained this way; these stellar masses are 
shown in Table [T] We make a final check of our method by 
comparing these derived stellar masses to those calculated by 
J. Huang et al., (in prep); this is shown in Fig. [T];, where the 
relationship is well-behaved with an rms scatter of 0.175 dex. 

2.8. Error Estimates 

All error estimates measured by GIM2D, i.e., the struc- 
tural parameters suc h as r^ and n, are 99% confidence limits 
jSimard et al.ll2002h ; we convert these into l-tr limits assum- 
ing a Gaussian distribution. There are two sources of stellar 
mass: those from the SED-fitting and those from our mass fits. 
The errors for the former are the typical widths of the stellar 
mass probability distribution (0.1 -0.2 dex). The errors for 
masses obtained from Eqn.|6]|3] and|4]are the standard devia- 
tion of the residual distribution between the fitted masses and 
those of J. Huang et al., (in prep) (see Fig.[T];; a w 0.175). Er- 
rors for U-B and Mb obtained from ^-correct are estimated 
by measuring the l-cr dispersion from the H ST /ACS V and 
/ in the redshift range 0.64 < z < 0.68 and 0.82 < z < 0.86, 
respectively. Within these redshift ranges, rest-frame U and 
B approximately redshifts into observed V and /, which when 
combined with the high resolution of HST, gives us an accu- 
rate photometric error estimate. The average errors are « 0.07 
mag (U) and « 0.05 mag (B). Errors for U-B and Mb ob- 
tained from Eqn.|3]and|5]are taken to be the standard deviation 
between our fits and ^-correct. For the mass-radius combina- 
tions, e.g., M,/re, we propagate the errors from the masses 
and the GIM2D confidence limits. 

3. SAMPLE SELECTION 

Within the AEGIS region, « 30,000 objects have both 
HST /ACS imaging and GIM2D decompositions; this is the 
master GIM2D sample. Only 11,223, however, have ei- 
ther a spectroscopic or photometric redshift (see More- 
over, although our redshift coverage is from 0.2 < z < 1.2, 
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Figure 2. Mf/Lg ratio (mass from J. Huang et al., in prep and B-band lu- 
minosity from i-correct) vs. observed V - I (left column) and rest frame 
U-B (right column) in three redshift bins: a)-b) 0.20 < z < 0.50, c)-d) 
0.50 < z < 0.80, and e)-f) 0.80 <z< 1.10. Contours are shown to give a 
sense of the relative number densities. The red dashed curve in each panel 
represents our fit; the fit for M* /Lg as a function of observed V — / is better 
than the fit to rest-frame U-B, hence we adopt it. 
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Figure 3. To address dust, we plot our sample in the rest-frame UVJ dia- 
gram. The quie scent population lie within the quiescent region as defined by 
IWilliams et aU 1 2009). Red points represent red sequence galaxies, defined 
to have U — B > 1.20. Almost all of the red sequence lies within the quies- 
cent area. U-B > 1.20 galaxies outside the quiescent area is only 3% of 
the defined red sequence. These are the dusty star-formers, and we eliminate 
them from the sample so that the red sequence galaxies are truly quiescent. 

in order to minimize ^-corrections, we restrict our sample to 
0.5 < z < 0.8; this cuts our sample down to 3,426, this will be 
referred to as the 'starting' sample. To reduce the effects of 
dust, we only choose galaxies with axis ratios {b/a; as mea- 
sured from the single n fit) greater than 0.55, furthering re- 
ducing our sample to 1,567 galaxies. 

Although GIM2D was run for every galaxy, not every de- 
composition is reliable. For example, galaxies with effec- 
tive radii r^ less than half the full width at half-maximum 
(FWHM) of the point-spread function (PSF; 2 pixels) are not 
well fit. Additionally, galaxy models created by GIM2D that 
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Figure 4. Le/f column : Rest-frame U—B vs. log M* is plotted for the 'starting' sample (red; top row) and the 'final' sample (blue; bottom row). For comparison, 
the Spitzer/IRAC 3.6 /im-selected sample (J. Huang et al., in prep) in our redshift range 0.5 < z < 0.8 is plotted in the background. Right column : U -B vs. 
log Mr is binned with lengths in U—B and log M, that correspond roughly to their l—cr error; we only show bins with more than 5 galaxies. Within each bin, 
the fraction of the number of galaxies in that row's sample to that of the 3.6 /im-selected sample is computed and displayed as the corresponding color indicated 
by the color bar to the right. The completeness of the 'starting' sample is uniformly complete above log M^/Mq = 9.5. The completeness of the 'final' sample 
is ~ 50% and is largely uniform. See text for discussion. 



are offset from the center of the HST/ ACS image by more 
than 3.5 pixels are similarly ill-fit. There are also instances 
where the fitting failed; eliminating these leaves us with 1,427 
objects. Note that we are only using the single Sersic fit values 
for global galaxy parameters, and hence this sample consists 
of values only from the single n fit. 

GIM2D was also used to produce measurements of every 
galaxy's bulge and disk through two different fits - the rthuige = 
4 and = 2 decompositions with the disk being n = 1 for 
both (see 32.31 ). Note that GIM2D bulge-i-disk decompositions 
do allow for a galaxy to have B/T = 0, i.e., a pure disk galaxy, 
if that is the optimal fit according to the Metropolis fitting 
alg orithm 22% of th e subcomponent sample have B/T = 0; 
see iSimardet al.ll2002l for more details). For each galaxy, we 
use the bulge-i-disk fit with the smallest x^- We only use the 
subcomponent measurements of the 'final' sample, which we 
define below. 

To reduce the effects of dust, we applied an axis ratio cut of 
b/a > 0.5 5.. While th is cut eliminates many edge-on dusty 
galaxies ( Martin et al.l 120 07). it does not affect dusty face- 
on galaxies. To clean these from out sample, we calculate 
t/y 7 rest-frame magnitudes. The resultant fZ-V vs. V-J 
two color plot enable us to separate dusty red galaxies from 
truly quiescent red galaxies (Williams et al. 2009). We use the 
Rainbow software described in Bar ro et alj(i2011a b). Briefly, 
the software applies a minimization algorithm to find the 
best fitting galaxy template from the multi-wavelength pho- 
tometry of AEGIS. Then several filters (U Bessel, V Bessel, 
J Johnson) are convolved with the best template to estimate 



synthetic fluxes assuming a luminosity distance of 10 pc. Our 
results can be seen in the UVJ diagram (Fig. O. The upper- 
left region bounded by the solid lines within th e UVJ diagram 
represents the quiescent region, as defined by I Williams et alj 
(2009). Comparing the quiescent galaxies to the red sequence 
galaxies, which we define to be galaxies with U-B > 1.20 
and are shown in red in Fig. [51 shows excellent agreement; 
only 17 (6%) of the U-B > 1.20 galaxies lie outside the 
quiescent area. These are presumed to be dusty, star-formers 
and are discarded from our sample. There is an additional re- 
duction of 8 fZ-B > 1 .20 galaxies because we were unable to 
obtain their UVJ magnitudes. Since we do not know whether 
these galaxies are truly quiescent or simply dusty, we take the 
conservative route and discard them. To sum, we require our 
quiescent galaxies to have U-B > 1.20 and to lie within the 
quiescent region of the UVJ diagram. With this criterion, our 
galaxy sample has 1,402 galaxies. 

3.1. Completeness 

Finally, we must discuss our sample's completeness. Be- 
cause the DEEP2-I-3 spectroscopic survey is limited by an R- 
band magnitude of 24.l[3, there is a selection bias against 
low-mass galaxies. Fortunately, the photometric sample goes 
deeper, down to an IRAC 3.6 /xm flux of 2 fi}y. Details 
of the photometric sample can be found in J. Huang et al., 

4 of these galaxies are not visible because they have UVJ magnitudes 
that are identical to those that are visible. 

There are some DEEP3 targets fainter than this limit ICooper et al.l 
I20nil20ll) 
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(in prep), but we will briefly summarize the key characteris- 
tics. The 3.6 /im-selected sample spans the redshift range of 
0.4 < z < 1.2, where 3.6 /im also probes the rest-frame NIR 
(1.2-2.5 fj,m). Galaxies of all types have very similar SEDs 
in the NIR band. Therefore a rest-frame NIR-selected sample 
suffers no b ias against ei ther blue or red galaxies ( Cowie et aTj 
[l996; Hua ng et al.lll997h . Galaxy NIR luminosities also trace 
their underlying stellar mass, in other words, this sample is 
very close to a mass-selected sample. The /T-band absolute 
magnitudes for galaxies in this sample are calculated with the 
3.6 /im flux densities. The IRAC-to-TT-band A:-correction is 
adopted from iDe Propris efaH dlOOTh . The absolute /T-band 
magnitude range for this sample is -19 < Mk < -25. This 
translates into a limiting stellar mass of log M^JMq =9.5. 
Therefore, cross-matching to the photometric sample has es- 
sentially eliminated the selection bias against low-mass galax- 
ies of the DEEP2H-3 surveys. 

To illustrate our sample completeness, we compare the 
color-mass diagrams of our 'starting' sample (red) to the 
SpitzerllRAC 3.6 /im-selected sample (gray) in Fig. |4^. 
There are hardly any gray points, indicating that the 'starting' 
sample contains almost all the galaxies in the 3.6 /i-selected 
sample. This is further illustrated in Fig.|^ where we bin up 
the color-mass diagram with lengths inU-B and log that 
roughly correspond to their distributions' l-cr error; we only 
show bins with more than 5 galaxies. Within each bin, the 
fraction of the number of galaxies in the 'starting' sample to 
that of the 3.6 /im-selected sample is computed and displayed 
as the corresponding color indicated by the color bar to the 
right. Confirming what was seen in Fig. 2^, the completeness 
is almost perfect, and most importantly, the completeness is 
uniform, especially for log M^/M© > 9.5, the mass limit of 
the 3.6 /im-selected sample. Thus, our 'starting' sample is 
uniformly complete down to the mass limit of the 3.6 /im- 
selected sample. 

However, our 'starting' sample is not the ultimate sample 
we use. To get rid of bad data and dusty galaxies, we have 
imposed several requirements (see To obtain our 'final' 
sample, we impose one final requirement, log M^/M© > 9.5. 
This last cut ensures that our 'final' sample is complete above 
log M^JMq = 9.5. Thus finally, we have our 'final' sam- 
ple, consisting of 943 galaxies. The 'final' sample is what 
is plotted in all subsequent figures unless stated otherwise. 
The completeness of the 'final' sample is illustrated in Fig.|4}; 
and Fig. |4}l. Fig |4}l (calculates bins of completeness like in 
Fig. lit) shows that the completeness of the 'final' sample 
is ^ 50%, with a dearth of galaxies on the top of the blue 
cloud, i.e., the green valley, and a surplus of galaxies on the 
upper red sequence. These features are due to the b/a cri- 
terion, which is meant to eliminate edg e-on galaxies t hat ar e 
presumably dusty. Indeed, according to iMartin et al.l (120071) . 
dusty galaxies do primarily reside on top of the blue cloud, 
which explains why there is a lack of galaxies on top of the 
blue cloud in the 'final' sample compared to the 3.6 /im sam- 
ple. The surplus of galaxies on top of the red sequence is also 
understandable since the reddest galaxies are elliptical galax- 
ies that have intrinsically high axis ratios. Although there are 
some biases introduced into the 'final' sample by these var- 
ious cuts, we have tested the effects of removing them and 
find that it does not affect our conclusions. But we stress that 
these cuts are necessary; they remove bad data. Our 'final' 
sample is a culmination of the best data from our available 
resources. For an extra discussion of our samples' surface 
brightness limits, data quality, and possible Sersic index bias. 



please see Appendix lAl iBl ICl 

All our data, including those that were not 
presented in this paper, are available online at: 
|htt p : / /people . ucsc . edu/~echeungl/data ■ html| 
Tables [l]|5] present the key parameters we use in our paper 
for twenty randomly selected galaxies in our catalog. Table 
[T] presents basic information of our galaxies, including 
their unique IDs, derived photometric quantities, and stellar 
masses both from A:-correct and Eqn. |6l |3] and |4] Table |2] 
presents much of the same information in Table 1, but only 
for the subcomponents. Tables |3] ID and |5] present the three 
GIM2D catalogs: the single Sersic fit, Ubuige = 4 fit, and 
rihuige = 2 fit, respectively. These GIM2D catalogs provides 
many measurements, including V, I, and for both the 
galaxy and its subcomponents. The rest of the measurements 
and galaxies can be obtained online. 

4. RESULTS 

4. 1 . The Most Discriminating Color Parameter 

We begin by intercomparing the various global structural 
parameters discussed in the introduction to find which is the 
best predictor of color. Fig. |5] plots U -B rest-frame color 
against six quantities for the final galaxy sample: rest-frame 
absolute B-band magnitude Mb, stellar mass M*, stellar mass 
divided by semimajor axis effective radius M^:/r^ (sometimes 
called the "inferred velocity dispersion"0 M, / (nominal 
surface densityS Sersic index «, and inner stellar m ass sur - 
face density Sij^p^ (we defer discussion of S*j,p^ to ^4.3.3b . 
The 1 -g error bars are given in the top right of each panel (see 
^2.8l for details). The spectroscopic sample and photometric 
sample are shown in open stars and open circles, respectively; 
we use this scheme throughout the rest of the paper As stated 
in the introduction, the amount of color overlap is one mea- 
sure of how well a parameter separates red sequence and blue 
cloud galaxies, and parameters that reduce this overlap are 
better discriminators of galaxy quenching. The sample con- 
sidered is the 'final' sample, which is complete only in stellar 
mass, as defined in ^ Thus the results of the following anal- 
ysis is only applicable for the 'final' sample. The goal of this 
section is to quantify the amount of overlap in order to deter- 
mine the single best color discriminant among the traditional 
parameters. 

As in IStrateva et al.l (|2Q()1|), the color-magnitude diagram 
(Fig. 5a) shows a clear red sequence and blue cloud. How- 
ever, these two groups overlap greatly over the entire range of 
absolute magnitude. This confirms the well-known result that 
the B-band magnitude is a poor predictor of galaxy color 

Fig. [5t) shows the color- mass diagra m. As shown by 
iKauffniann et all (l2003h and lBorch et all (l2006h . mass is bet- 
ter correlated with star-formation history than is luminosity. 
Although the relationship with color is improved, the range of 
color overlap is still large, extending over ~ 0.8 dex in mass. 
Fig. |5}; and|5}l add powers of r^ in the denominator to M,, in 
the form of M^/r^. The smallest overlap by eye is given by 

The true stellai' velocity dispersion is oc GM/r^, wliere M is tlie to- 
tal mass including stars, gas, and dark matter Franx et al. ( 2008) provide 
a value of the coefficient through the fitting of a sample of SDSS galaxies: 
0-2 = 0.3GM,/re. Recendy, Taylor et al. (2010) and Bezanson et al. (201l|) 
showed that the addition of a Sersic dependent term to the "inferred velocity 
dispersion" of Franx et al. (2008) provides a better estimate of the true veloc- 
ity dispersion. We choose not to use this updated "inferred velocity disper- 
sion" because we want to compare the color correlations of these parameters 
independently. 

Surface mass density is actually Mt/lnrl, but we omit the constants. 
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Figure 5. Rest-frame U —B color is plotted against: a) rest-frame absolute B-band magnitude Mb, b) stellar' mass M*, c) stellar mass divided by semimajor axis 
effective radius M-t./r^, d) M^/rl (surface density), e) Sersic index n, and f) stellar mass surface density within 1 kpc 2*|,pj.. M* is in solar masses and re is in 
kpc. The stars and circles represent galaxies with spectroscopic redshifts and photometric redshifts, respectively, l-a error bars are displayed in the upper right 
of each panel. Th e ove rlap region of each parameter is highlighted in gray. Overlap regions are defined as regions in which the fraction of red galaxies is between 
15% -60% (see j|4. It . The red horizontal line represents the division between red and blue galaxies. The top left comer of each panel shows two measures of 
the size of the overlap region: the fractional number of galaxies within the overlap region AN/N and the fractional extent of the region Ax/x, wh ere x contains 
90% of the points (see text). This analysis was done on the 'final' sample, which is complete only in stellar mass (down to log M» /Mq = 9.5; see j|3.1t . By both 
measures. A/* /r^ outperforms other parameters of the form A/* /r^ , which we have confirmed by studying intermediate values of p. Sersic index and S^j^p^ do 
even better, hinting that the distribution of mass in the inner parts of galaxies may play a fundamental role in quenching star formation. 



M*/re in Fig.|5};, while mass surface density M^,/rl in Fig.|5}l 
looks slightly worse. Thus, this new DEEP2 sample indicates 
that M^,/rg is a superior color discriminator to surface mass 
density M^,/rl. 

To summarize, effective radius tightens the basic color- 
mass relation because red galaxies at fixed mass are smaller 
than blue galaxies, and the tightest correlation is obtained us- 

We plot in Fig. color vs. Sersic index n. The charac- 
ter of this plot is markedly different from the others - rather 
than a smooth trend with color within the blue cloud as in, 
for example, M^jr^, the color jump is more abrupt, with 
color remaining constant above and below what appears to 
be a critical value of n around log n = 0.36 (« = 2.3). This 
behavior is intriguing because it might signal a real phys- 
ical threshold in Sersic index, above which star formation 
shuts down. As stated in the introduction, n likely plays 
an imp orta nt role in quenching star-formation. Blanto net alJ 
(|2003) and ISchiminovich et alJ (120071) demonstrated a trend 
bet ween n and color for SDSS galaxies. Driver et al. (200^ 
and lAllen et al.l (2006) also showed this r elation ship with their 
Millennium catalog. And IWuvts et"an (1201 1[) demonst rated 
that this relation ship persists out to z ~ 2.5. iBelll (120081) and 
IB ell et af] (120121) explored this correlation and showed that 
high n is necessary for quenching but not sufficient - there 
are many galaxies that are blue despite having high n. We 
see something similar in our data with the scattering of aber- 
rant points in the lower-right-hand comer of Fig. |5^. We term 



t hese a berrant points "outliers" and discuss them further in 
ESS 

The above conclusions are based mainly on visual assess- 
ment of Fig. |5] To quantify our results, we now present two 
new quantities that are designed to measure the size of the 
overlap regions. These measurements can be applied to rank 
the predictive power of the various structural parameters and 
also to identify galaxies within the overlap regions for fur- 
ther study. To define these quantities, we first bin the sam- 
ple by the parameter of interest. Then, within each bin, we 
find the fraction of galaxies that are red, i.e., galaxies with 
U-B > 1.20, which we have ensured to be genuinely quies- 
cent and not dusty (see The locations where this fraction 
equals 15% and 60% mark the beginning and end of the over- 
lap region, respectively; these percentages were adjusted to 
match the core of the overlap regions as judged by eye and 
are a compromise over all diagrams. We varied the overlap 
definition with various permutations of starting boundaries in 
the red fraction range of 5% to 20% and ending boundaries 
from 40% to 60% and found that the results are unchanged. 
To find the locations of the edges of the overlap regions, we 
fit a fourth-order polynomial to the red fraction bins and in- 
terpolate to find where the fit reaches the desired fractions. 
Each parameter has been divided into 25 bins, and each edge 
value is examined to ensure that the choice is sensible. The 
edge locations depend only weakly on the choice of bin width 
- for example, in the case of M^,/r^, bin sizes in the range 
0.07-0.25 dex produce similar results. 
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We define two measures to quantify tlie sizes of the over- 
lap regions. The first is the fractional number of galaxies in 
the region, AN/N, where N is the total number of galaxies and 
AN is the number within the overlap region. The second is the 
fractional extent of the region, Ax/x, where Ax is the width 
of the overlap and x is the range that includes 90% of the 
data (excluding 5% at either end). The resulting overlap re- 
gions for each parameter are demarcated in gray in Fig.|5] and 
the upper left corner of each panel shows the two measures 
AN/N and Ax/x. These quantitative measures confirm what 
was seen by eye, namely, that M^jr^ gives the smallest values 
of both AN/N and Ax/x among all the mass-radius combina- 
tions. Note that the relative extent of Mb is Ax/x=0.971; this 
simply means that the overlap region is almost equal to the 
entire range that encompasses 90% of the data, again agree- 
ing with our previous qualitative assessment. Furthermore, 
we find that Sersic index performs considerably better than 
even M^jr^ in minimizing both AN/N and Ax/x. 

We also point out the extremely tight relation that is pro- 
duced when plotting color vs. M^jr^ or for blue-cloud 
galaxies alone (Fig.E}? and|5]:). This has been pointed out be- 
fore and is referred to as the 'Main Sequence' of star forma- 
tion (jNoeske et al. 2007). Previous work on quenching has fo- 
cused on the relationship between red sequence vs. blue cloud 
galaxies and not so much on the properties of galaxies within 
the blue cloud itself. However, the tightness of the relation 
between M*/re (or M*) and color for star-forming galaxies 
alone could be an important clue to the physics of quenching, 
and we return to this point in ^ 

4.2. Properties of Galaxies in the M^/r^ Overlap Region 

We have shown that our AEGIS data duplicate previous 
findings showing that M^Jr^ and n correlate strongly with 
quenching, but we have also shown that neither parameter 
alone is close to being a perfect predictor of it. In this sec- 
tion and the next, we take a further look at the properties of 
galaxies in the overlap region and outliers to find out whether 
multiple parameters can be used in concert to predict quench- 
ing, and whether this interplay sheds light on the physical pro- 
cesses involved. 

Fig.|6]investigates residual trends within the M^jr^ overlap 
region of Fig. |5}; by plotting color versus various structural 
parameters for overlap galaxies alone. In exploring this slice 
of M*/re, we are assuming a general evolutionary sequence 
such that the blue galaxies evolve into the red galaxies within 
this overlap region. However, this assumption is not without 
proof. Bell et al. (2004b) and Faberet al. (2007) have shown 
that the red sequence has increased by ~ 2 while the blue 
clou d has remained relatively stable from z 1 to z = 0. More- 
over, iHopkins et al] (l2010bh showed that 65 - 80% of the 
observed mass density of bulge-dominated galaxies formed 
since z ^ 1. Thus our redshift range (0.5 < z < 0.8) peers di- 
rectly into the epoch when the majority of red sequence galax- 
ies are being formed. We choose M, / r^ as the base parameter 
because it is the tightest M, - r^ combination in Fig. |5] but 
similar results are obtained when M, is used. 

Panels a, d, j, and k of Fig. |6] plot integrated quantities, 
while the remaining panels introduce structural parameters 
(e.g., M, and r^) for bulges and disks separately. Among the 
integrated properties, virtually no trend is seen in either stel- 
lar mass (panel a) or r^ (panel d), but a mild "step function" 
is seen with Sersic index in that n is significantly higher for 
quenched galaxies (panel k, gray). Concentration, C, shows 
similar behavior, albeit less cleanly (panel j). A similar trend 



with n was seen for all galaxies (Fig. |5^), but it is important 
to see the same effect for overlap galaxies alone. This estab- 
lishes beyond doubt that M^jr^ alone does not fully encapsu- 
late the processes needed to quench star formation. A possible 
interpretation is that some galaxies are "ripe" for quenching 
based on M, / r^ and that a second process, which drives galax- 
ies to high n, ultimately quenches star formation. We return 
to this idea later in ^ 

The remaining panels in Fig. |6] focus on the properties 
of bulge and disk components separately (the subcomponent 
sample; see ^2. 3b . These parameters are derived from G1M2D 
photometric fits and M,/Lb values from the Y -I colors of 
bulges and disks separately 02.7l i. for which high-resolution 
two-color HST imaging is required. The striking result from 
the subcomponent panels in Fig. |6]is the marked differences 
in disk mass, bulge mass, and bulge effective radius between 
blue and red overlap galaxies (panels b, c, and e). The disks 
of red sequence galaxies are less massive by about 0.2 dex 
than the blue cloud galaxies (panel c) while the bulges of red 
sequence galaxies are more massive by about 0.3 dex than the 
blue cloud galaxies (panel b, in light tan). At the same time, 
disk sizes remain constant but bulge sizes decrease by about 
0.3 dex as blue cloud galaxies transition to the red sequence 
(panel e, in blue). These differences between red sequence 
bulges and blue cloud bulges amount to a difference of 0.6 
dex in bulge bulge /'"e.buige (panel h, in green). A similar 
trend is seen m B/T ratios, but it is weaker due to the large 
spread of blue galaxies. 

These trends collectively demonstrate a real structural 
difference between the inner stellar mass distributions of 
quenched vs. star-forming galaxies, and furthermore, that this 
difference exists even within a narrow range of / r^ . Higher 
central stellar mass densities in quenched galaxies have been 
inferred in previous work from higher integrated Sersic val- 
ues ( Weiner et al. 2005; Bundy et al. 2010), but photometric 
data by themselves do not rule out a simple fading picture in 
which disks decline in brightnes s, permitting an un derlying 
high-« bulge to emerge (see also iHolden et al.ll2012l for dis- 
cussion of b/a). Actual stellar masses for bulges and disks 
separately are needed to rule out fading. An important new 
insight from our work is that evolution to the red sequence 
appears to be accompanied by a significant rearrangement of 
inner stellar mass in which existing stars move to the centers 
of galaxies, and/or new stars are formed there. We discuss 
processes whereby that might happen in ^ 

4.3. Sersic Index and Inner Surface Density 
4.3.1. A threshold in n? 

The previous section considered M, / r^ as the main quench- 
ing parameter and looked at scatter around it to discover sec- 
ondary effects. In this section, we use a similar tack but focus 
on Sersic index n. An enlarged version of the color-« plot is 
shown in Fig.|7^, which indicates median U -B in bins of n. 
The medians illustrate the striking step-like behavior previ- 
ously mentioned. Defining the riser of the step to be where 
the medians are half-way between their red and blue values, 
we place this point at log n = 0.36, which is marked with the 
vertical black line. This value corresponds to n = 2.3, which is 
similar to the value of n often used to distinguished quenched 
(or early-type) galaxi es from star-forming ones locally (e.g. , 
Shen et al. 2003; Bell et al.ll2004at ISchiminovich et al.ll2007l: 
iDrorv & Fisheiii2007i) . 

The medians also emphasize how flat the color trends are 
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Figure 6. Galaxies in the overlap region of Mf/r^ (containing 316 galaxies) from Fig.[5j; are isolated to examine possible second-parameter color coiTclations 
with both global and disk/bulge properties. U —B rest-frame color is plotted against: a) global stellar mass M, , b) bulge stellar mass M« bulge, c) disk stellar mass 
^•,disk, d) global semimajor axis effective radius r^, e) bulge semimajor axis effective radius J'cbuigc. f) disk semimjaor axis exponential scale length jis^, g) 
^*,bulge/''e bulge' bulge /''e.bulge, bulge-to-total ratio B/T in the /-band, j) concentration C, k) global Sersic index n, and 1) B/T based on stellar mass M*. 

The l-CT error bars are displayed in the upper right of each panel. The red horizontal Hne in each panel represents the division between red and blue galaxies. 
Global properties hardly differ between red and blue galaxies (white panels), but bulge properties, disk scale length, and n show key differences (colored panels 
and panel c). Specifically, red sequence bulges are physically smaller, yet more massive, than those of the blue cloud. Moreover, red sequence disks are less 
massive than their blue cloud counterparts. These trends cannot be produced by simple disk fading but require a concentration of inner stellar mass. 

low and high redshifts by iBelll (120081) and iBell et af] (120121) . 
who expressed the role of n in quenching as "necessary but 
not sufficient", i.e., all quenched galaxies have high n, but not 



above and below the threshold n value. Evidently, in the ex- 
treme high-n and low-« regimes, star formation history is in- 
dependent of n. This differs markedly from the behavior of 
M^,/r^\ Fig.|5}; shows a strong correlation between /rg and 
color for star-forming galaxies below the overlap region. 

The lack of importance of n above and below the thresh- 
old is further emphasized by the large color scatter in both 
of these regimes. This scatter is, however, of two types. At 
low n, there is a rather uniform spread in color, i.e., specific 
star-formation rate can assume any value within a large range, 
and n does not predict what SSFR will be. At high «, n pre- 
dicts color much better, i.e., the color distribution is strongly 
peaked toward red (quenched) colors, but a significant tail of 
outliers with blue colors exists (colored points, except the red 
in Fig. |7ll. The existence of these outliers was seen at both 



all high-« galaxies are quenched. We see the same thing. 

Unlike Fig.jS}; (which plotted color vs. M„/re), there is no 
interval in n where the color scatter is markedly larger than 
elsewhere (Fig.|7^), and thus no impetus to search for a sec- 
ond parameter within a narrow region of n. To investigate 
the scatter, we have replotted Fig. |7] this time highlighting 
galaxies within narrow bins of various second parameters. 
The results are shown in Fig. [Sj where galaxies are divided 
into bins of stellar mass (top row), M*/re (middle row), and 
M^,/rl (bottom row), collectively termed M^./r'e. The outlier 
region from Fig. is outlined in blue. In each row, the be- 
havior is the same. Galaxies with low values of /r^ are 
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Figure 7. a) The U — B vs. n diagram is reproduced from Fig.|5^ with median values of color in each n bin plotted as open black triangles. The red horizontal 
line represents the division between red and blue galaxies, while the black vertical solid line at log n = 0.36 represents the point at which the median is half way 
between the red and blue values. The medians show a step-like behavior in U —B vs. n, suggestive of a physical thres hold in n. Exceptions to the step function 
are the blue "outlier" galaxies in the lower right. These are highlighted in strips of color for further discussion in j|4.3.2l and Fig.|9] pink points encompass 
U —B < 0.70, cyan points encompass 0.70 < U -B < 0.95, and green points encompass 0.95 < U -B < 1.20. Bluish red-sequence galaxies ai'e highlighted for 
comparison and lie within 1.20 < U —B < 1.30. Images of these outliers are shown in Fig.|9] Roughly 40% of the log n > 0.36 galaxies are outliers, b) U — B 
color vs. inner stellar mass surface density Note that some galaxies do not have S*j^p^ measurements due to insufficient signal-to-noise. Most outliers now 

fall into line, suggesting that inner stellar mass surface density is a cleaner predictor of quenching than n. 



seen to be mainly blue. A few leak into the high-« "outlier" 
regime, but their blue colors always agree with other galax- 
ies in the same parameter bin, i.e., their star formation rates 
are not depressed by having high n. As M^,/f{ increases, the 
mean color of low-« galaxies becomes redder while the num- 
ber of outliers remains relatively constant. Again, the colors 
of the outliers still agree with the average color of all galaxies 
in the same M^,/r'l bin. At the highest values of M^ jr'i, vir- 
tually all galaxies are quenched and the fraction of outliers is 
negligible. Two points are clear: dividing galaxies into bins 
has not destroyed the basic step-like nature of the behavior in 
that galaxies within each bin still trend smoothly but sharply 
(apart from outliers) from their "native" star-forming state to 
a quenched state. The second point is that all trends with color 
at low n remain flat within each bin of / r'l . This shows that 
the strong trend of color vs. or / within the blue cloud 
(Fig El' and Eh) is not caused by some hidden dependence on 
n. 

To summarize, we have reproduced findings by previous 
authors that indicate that high n typically predicts a quenched 
galaxy, and we have set the half-power point between blue and 
red galaxies at = 2.3. This value is very near the SDSS 
value, implying no large evolution in from z ^ 0.65 down 
to z ~ 0. The rise in color near the critical value is sharp, while 
above and below this value there is no trend in color with n, 
even within narrow slices of M,/re . At high n, most galaxies 
are quenched with red colors, but a non-negligible fraction of 
objects is blue despite having high n. We turn to the nature of 
these outhers next. 

4.3.2. Outliers 

Although n acts like a threshold for the vast majority of 
galaxies, there are obvious exceptions, namely the blue, high- 
n (logn > 0.36) "outliers" highlighted in Fig. and else- 
where. Understanding these objects is clearly crucial for un- 
raveling the quenching mechanism - why are they blue when 
their photometric structure resembles that of quenched ob- 
jects? We have highlighted 151 outliers in Fig.|7]using color 
to indicate their U -B ranges; they make up w 40% of all 
n > 2.3 galaxies (the red points immediately above the red 



horizontal dashed line at U -B = 1 .20 are not outliers, they 
are quiescent red sequence galaxies shown for comparison). 

Several possibilities come to mind to explain these ob- 
jects. One possibility is that they are artifacts due to the 
presence of bright point-like AGNs. Adding a blue AGN 
to a normal star-forming galaxy would make th e global col- 
ors b luer and increase n (and concentration) (iPierce et al.l 
I20T0I) . To pursue this, we have cross-matched the outliers 
to two AGN samples selected using X-ray and optical line- 
emission data. The AEGIS region is covered by a deep 800 
ksec Chandra X-ray mosaic (iLaird et al.l (2009). We find 
that only 11 of the outliers have X-ray luminosities above 
10'*^ erg s"', or 7%. We have also used an optical method 
for selecting AGN based on a modified version of the classi- 
cal "BPT" diagram (Baldwin, Phillips & Terlevich 1981) that 
plots [O111I/H/3 versus U -B rather than [Nll]/Ha (lYan et all 
1201 II) . This adds only 14, bringing the total to 21 AGN, or 
14% of the outliers. Thus, it seems that the vast majority of 
these objects are unlikely to be AGN hosts. 

Another possibility is a bright blue clump of recently 
formed stars at the centers of the outliers, which might skew 
the colors and Sersic indices as an AGN would. We would like 
to remind the reader that the single Sersic G1M2D model does 
not fit for any substructure, such as clumps. Hence asymmet- 
ric structures may affect the n measurements. To explore this, 
we constructed V and / color images of all outliers using the 
HST/ ACS data and inspected them; a montage of 20 galaxies 
is shown in Fig.|9] where each row represents a different range 
ofU-B color according to the color-coding in Fig.|7^. The 
bluest outliers are in the bottom row. These tend to be lumpy, 
asymmetric, and/or small - their fitted Sersic values are some- 
what questionable. Moving up one row to the cyan points, 
we find a mix of genuinely concentrated galaxies plus more 
small galaxies like the ones in the previous row. The third row 
contains larger objects of intermediate color but with convinc- 
ingly high concentrations. Finally, we show a sample of red 
sequence galaxies in the top row as a comparison; they are all 
red and highly concentrated spheroids. 

To summarize, the blue, high-n outliers are a mix of dif- 
ferent types. Some may have doubtful Sersic indices, be- 
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Figure 8. The U — B vs. n diagram from Fig.|7^ is replotted to separate the roles of M«, Mtjr^, and M^/rl from n in driving color evolution. Galaxies are 
divided into bins of M, Ir'l and plotted as the black points. Bin boundaries are shown in the upper left of each panel. The rectangles outlined in blue indicate the 
outlier region from Fig.|7^. The panels demonstrate that the step-function-like behavior near log n = 0.36 (n = 2.3) seen in Fig.|7p is replicated separately in each 
range of M, /r^ independent of p. Color rises rapidly near log n = 0.36 (n = 2.3), but there is little impact of n on color above and below this value. This shows 
that the trends in color with M* and M» / in the blue cloud (Fig.|5j) and|5];) are not caused by a hidden dependence on n. 



ing small or with off-centered clumps of star formation or 
AGN, but a fair fraction seem to be genuinely blue yet high- 
n galaxies. These genuine exceptions tend to be located at 
intermediate values of U - B. The existence of such out- 
liers has been noticed before. An i mportant class of candi- 
dates is poststarburst ga l axies (e.g.. iDressler & G un nlll983l; 
Couch & SharplesI 119871: iPoggianti et al.1 11999: Go to et all 
2003"). These objects possess blue colors and strong Balmer 
absorption yet weak Ha, sign ifying recent quenching, and 
their Sersic indices are high (Ouintero et al. 2004^, lYang et all 
|2008). A second class of objects is the blue spheroidal galax- 
ies; like poststarbursts, they have smooth, centrally concen- 
trated, elliptical-like profiles but they are different in having 
active star forma tion (Menanteau et al. 2001; Im et al. 2001; 
Koo et all 120051: ISchawinski et al.i |2009; Kannappa n et all 
2009^ . Their masses tend to be small (Imet al. 200 C and 



there appear to be several objects in our outlier population 
that fit this description in the bottom row of Fig. |9] 

A quick calculation of the percentage of outliers within 
a volume-limited SDSS sample at z ~ shows that it has 
dropped from w 40% for our sample at z ^ 0.65 to « 10% 
at z ~ 0. This difference seems consistent with the higher lev- 
els of gas at higher redshifts, which could give rise to more 
clumpy star formation asymmetrically distributed throughout 
the galaxy, skewing the Sersic values. 

For completeness, in passing we also mention satel- 



lite processes. Pro cesses such as r am pressure stripp ing 
( iGunn & Gottll972l) and strangulation (iLarson et all 19801) do 
not by themselves change n. However, other satellite mecha- 
nisms involving tidal interactions (such as "harassment" e.g., 
iMoore et al.llT996 ) may induce a morphological transforma- 
tion. If satellite galaxies are first harassed, they might attain a 
high n while still forming stars. While we do not expect most 
of the DEEP2 galaxies to be satellites (see introduction), a 
more thorough investigation is needed to confirm this. 

4.3.3. Inner Surface Density: An Improvement on n? 

From the above, it is clear that in n we have found a re- 
markable, but still imperfect, structural predictor of quench- 
ing. The main criticism of n is the presence of outliers; if they 
could be removed, the correlation, and hence the prediction, 
would be nearly exact. 

There are two obvious weaknesses with n. First, it is based 
on light, not stellar mass, and hence is subject to the vagaries 
of star formation history and dust. Second, it relies on a fit to 
the entire light profile and is thus at least partially dependent 
on the outer light distribution, which may be disturbed or ir- 
regular In contrast, trends discovered using bulge properties 
in Fig.|6]hint that the structure of the inner regions of galaxies 
is more important at predicting quiescence. 

Accordingly, we introduce a new parameter that attempts 
to remove both of these weaknesses, namely the stellar mass 
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Figure 9. H ST /ACS V and / color images for a selection of high-H outlier galaxies from Fig. [TJi. Objects are an'anged in rows according to U —B color with 
bluer galaxies at the bottom. The bottom row corresponds to the pink points in Fig. [Tji with U -B < 0.70; the next row up coiTesponds to cyan points with 
0.70 < U-B < 0.95; and the following row is the green points with 0.95 < U-B < 1.20. The top row contains red sequence galaxies shown for comparison; 
they are represented by the red points in Fig.|7^ and have 1.20 < U-B < 1.30. A scale of 5 kpc at the average redshift (z si 0.68) of our sample is shown as a 
reference in the upper right picture. Many of these blue high-n objects seem to have Sersic indices perturbed by small size, irregular structure, central starbursts, 
and/or AGN and are cured as outliers if inner stellar surface density 2*|^pj, is used instead of n, as shown in Fig.|7j). 

GIM2D fit0 and that using a more secure quantity like in- 
ner mass surface density can remove them. It is still true that 
using Sii^pj. by itself is not perfect and that n somewhat out- 
performs it on the overlap criteria seen in Figs.|5^ and|5]F. That 
a genuine spread exists in S^'j^p^, is confirmed by ongoing work 
with higher-S/N SDSS data (Fang et al., in prep), which how- 
ever reveals some additional striking regularities. Our point 
for now is that using '^1^^^ removes the large number of out- 
liers that are present with n. Furthermore, the definition of 
^tkpc ™ inner mass density directs our attention even more 
strongly to the fact that it is conditions near the center of the 
galaxy that drive quenching. 

5. DISCUSSION 

In this paper, we have found that the most discriminating 
parameter of quiescence, according to the measures of the 

GIM2D only models a galaxy into either a bulge+disk model or a single 
Sersic model. More complex structures like spiral arms, bars, and clumps are 



surface density within 1 kpc. This is defined as Sji^p;, = 
■'^rkpc/'''''"ikpc' where rikpc 1 kpcQ The actual diameter used 
is 12 pixels (0.36"), which is set by the smallest radius that 
our HST images can conveniently resolve; it spans a radius 
of 1.08- 1.35 kpc at our redshifts. The quantities / and V -I 
are measured within this aperture, and Lb, M^,/Lb, andM* are 
estimated using Eqns.[3]|4l and|6l 

The quantity S*|jp(, was already included for completeness 
in Fig.|5](panel f), where its performance is seen to be mixed. 
It seems to predict color less well than the global quantities 
and M^,/r^ in the blue cloud, but it does much better 
than n in removing outliers. This is even better illustrated in 
Figs.|7^ and|7j5, which compare outliers directly. Using S*i5p(,, 
almost all the pink and cyan points have receded back into the 
blue cloud, and only a few green outliers remain. This sug- 
gests that the outliers in n are largely artifacts caused by poor 



^' Note that here we include the tt, so these are physical surface densities. not modeled. Thus a galaxy with these features are potentially ill-fit. 
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overlap region, is the Sersic index, and that a plot of color vs. 
Sersic index shows a step-like behavior near n = 2.3, sugges- 
tive of some sort of genuine quiescence threshold. About 40% 
of the n > 2.3 galaxies, however, are "outliers" that fall out- 
side this behavior These outliers have central mass densities 
much lower than those of the red sequence and fully consis- 
tent with those of the blue cloud. In other words, under this 
new parameter, Sji^p;,, the outliers now fall in line, suggesting 
that a galaxy's central structure may be even more physically 
related to quiescence than Sersic index. 

Both n and T^*^^^,^ corroborate our second major finding, that 
most blue cloud galaxies at the observed epoch cannot simply 
fade onto the red sequence. We have shown through stellar 
mass measurements of bulges and disks that the red sequence 
galaxies have absolutely higher bulge mass concentrations, 
i.e., t hat th e jump in n is not due simply to the fading of disks 
(see 34.2l i. The central mass densities extend this conclusion 
to the very centers of galaxies. In other words, the transition 
to the red sequence involves a significant restructuring of a 
blue cloud galaxy's innermost stellar density profile. 

Below the critical value of n = 2.3, however, Sersic index 
shows little to no correlation with star formation, and color is 
more closely correlated with M^/rg (or perhaps with M*, see 
Fig.E). 

This two-pronged behavior suggests that the star-formation 
history of a galaxy is shaped by two separate factors at 
different stages. While the object is still star-forming (in 
the blue cloud), its star-formation rate depends on global 
parameters, like M^,/rg or M^,. Then, a major internal mass 
reorganization occurs, a dominant bulge forms, and star 
formation stops. In the following discussion, we compare 
these results to the predictions of the standard merger 
model for bulge-building and quenching and find reason- 
able agreement, but also several issues. To alleviate these 
issues, we also consider other models, specific ally, violent 
disk instabihty (iNoguchil Il999t Elmegreen e tall 120081: 
Dekel et ani2009 '). secular evolution dKormendv & Kennicutt 
2004ft . morphological qu enchin g (Martig et al. 2009), 



and halo quen ching (Silk 1977; Rees & Ostriker 1977; 



Bhimenthal et al. 1984; Birnboim & Dekel 2003; Keres et al.. 
2005; Dekel & Birnboim 2006; Cattaneo et al. 2006). We 
end with a brief discussion on an unresolved concern. 

5.1. Merger Model 

Major merg erf-^ have been linked to the formation of 
spheroids since iToomre & Too mre ( 197 2|), with con siderable 
work in the years since (see IHopkins et al.l l2009al and ref- 
erences therein). The process of bulge formation in classi- 
cal merger models occurs through both the violent relaxation 
of pre-existing stars to the center a nd a gaseous dissipation- 
induced nuclear starburst (Hopkins et al. 2009a). 

Comparison of the bulge-dominated products of these 
simulations to observed early-type galaxies shows good 
agreement. The Sersic indices of the merger products 
from [H opkins et al. (2008) are high, n > 2.5. Similarly, 
iNaab & Truiillol (1200 6) showed that their disk mergers (with 
bulges) created galaxies with 3 < n < 4 [3- Both of these 



According to IHopkins et alj )2010at) , major mergers dominate tlie for- 
mation and assembly of ~ L« bulges and the total spheroid mass density. 
Thus, we only consider major mergers in this dis cussion. However, it sh ould 
be noted that minor mergers can create bulges <Boumaud et alj|2007l) and 
do contribute a non-negligible amount (~ 30%) to the total spheroid mass 
density fHookins et al. 2010a). 

iNaab & Trujillo ( 2006.) conducted a collisionless simulation that does 



works produce spheroids with Sersic indices in the range 
of ou r red sequence spheroids and of other observations 
(e.g., 'Kormendv & Kennicutt 120041; iDrorv & Fished 120071; 
Fisher & Drory 2008). 

The properties of quenching induced by mergers are 
also consistent with our data. The merger model pre- 
dicts that quenching oc c urs t hrough the nuclear star- 
burst dMihos & Hernauis3 119941) . in which a large por- 



tion of t he gaj^ is exhausted (depending on gas frac 



tion; see Soringel & Hernauist 2005; Okamoto et al.' '2005t 
Robertson et al. 2006; Cox et al. 2008; Governato et al. 200j, 
120091) . and throug h subsequent feedbac k, from super- 
novae and/or AG N dSpringel et anr2005; Murray et alj|2()0l 
ICiotti et £11120091) . Note that both sources of quenching orig- 
inate from the center of the galaxy, suggesting that the con- 
ditions at the center may correlate better with the quench- 
ing state than global properties. This is what we find. Fur- 
thermore, the central surface mass densities of the simulated 
spheroids from lHopkins et af] (12009b ) match our observations 
of the red sequence spheroids in Fig. [TJ) - with values of 
log S = 9 - 10 M0kpc-2 at 1 kpc. 

A further success of the merger hypothesis is the good 
match between it and the stellar mass range where bulge- 
building is observed to occur. A major point is that the ef- 
ficiency of bulge-building from major mergers is expected 
to be highly dependent on the pre-merger gas fraction, such 
that decreasing gas content increases the potential to form 
bulges (Springel & Hernquist 2005; Robertson et al. 20061 
IHopkins et al.ll20 09a). This dependence is consistent with the 
assumption that gas content gradually falls as galaxies age in 
the blue cloud, making them ultimately ripe for spheroid for- 
mation via mergers. Because more massive galaxies exhaust 
their gas qui cker due to the p henomenon known as "down- 
sizing" (e.g.. lCowie et alj |1996). there is a strong color-mass 
relation in the blue cloud, meaning that the reddest blue cloud 
galaxies at any epoch have the least amount of gas. Accord- 
ing to the merger model, this means that they are also on the 
threshold of forming bulges. 

Evidence for this hypothesis is strong in Fig. ^p, which 
shows a remarkably tight correlation between stellar mass and 
color in the blue cloud in the sense that more massive galax- 
ies are the reddest. Further data come from ICatinella et 
(120101) and Saintong e et all (1201 Ih . who presented H I and 
CO data in the GASS and COLD GASS survey, respectively. 
These works illustrate that the average atomic and molecu- 
lar gas fraction of galaxies do decrease with increasing stel- 
lar mass and increasing NUV-r color Although these sur- 
veys do not fully sample the blue cloud (these surveys only 
observe > 10'" M© galaxies), extrapolating these seem- 
ingly linear trends to lower masses indicate that total gas frac- 
tion does indeed decrease with mass along the blue cloud. 
Theoretically, Hopkins et al. (2010a) showed that the most 
effective bulge-building major mergers are clustered around 
log ~ 10.5 Mq at z ~ 1, which is in the center of over- 
lap region of M*. This behavior is due in part to exactly the 
same reason, namely, decreasing gas content as galaxies age 
within the blue cloud. In a general way, then, theory predicts 
that galaxy colors and gas contents should both age within the 



not include gas, and thus does not model a nuclear starburst component. That 
is why their pure disk-disk mergers only have « 1.5, because they lack the 
high central surface brightness typical of dissipative gas-rich mergers. 



35 



In the following discussion, gas is assumed to be cold gas unless other- 



wise stated 
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blue cloud, causing galaxies to become more prone to bulge- 
building mergers at higher mass and low gas level, and these 
trends broadly agree with the observations. 

In detail, however, the data indicate that M^,/r^ is a better 
predictor ofguenching than M, alone (cf. overlap regions in 
Fig. I5J3 vs. ISp). T his may be becau se M^/r^ is related to ve- 
locity dispersion jFranx et al.ll2008t see footn ote 26), which, 
based on a new studv bv IWake et al] (l2012bl) . is the galaxy 
property most closely related to halo mass. This finding could 
then be a manifestation of the dependence of quenching on a 
critical halo mass. Alternatively, it may reflect the fact that 
radii are shrinking as stellar mass builds up in the centers of 
quenching galaxies and thus reflects a property of the galaxies 
themselves rather than of their halos. We elaborate further on 
these thoughts in ^5.51 

Although there are some aspects of the merger model 
that match our data, general agreement upon the valid- 
ity of this model has not been reached. One controver- 
sial issue is whether major mergers can actually quench 
galaxies. There have been various wor ks that support 
this idea (e.g., [S chawinski et al. 2007; Alexa nder et al.ll2010l : 
Cano-Diaz et al. 2011; Farra h et al.i i201 ll) . In particular, 
Cano-Diaz et al. (2011) obtained VLT-SINFONI integral- 
field spectroscopy for one quasar at z = 2.4 and showed a sup- 
pression of narrow Ha emission, a tracer of star formation, in 
the region with the highest outflow velocity and highest ve- 
locity dispersion. But this is only one example and does not 
erase the contradicting evidence others have offered. For ex- 
ample, using a sample of X-ray and post-starburst galaxies 
from SDSS and DEEP2 at 0.2 < z < 0.8, C oil et al. (2011.) 
found winds with velocities that are insufficient to shut down 
star formation, indicating that the presence of an AGN does 
not produce faster winds nor does i t seem to play a major role 
in quiescence. And' Ammons et al] (1201 1 ) fail to find any cor- 
relation between host galaxy color and X-ray hardness ratio 
among z = 0.5- 1.5 galaxies, as might be expected from the 
blowout model. 

Another important concern is whether there are enough 
mergers to account for the bulge density in the universe. 
Studies that have measured the galaxy merger rate often 
present different results (e.g .. Bell et al. 2006; Conselice 
2006t iLin et a .Il2008t iLotzet ak 2008b.; Bundy et al. 2009; 
Xu et alJ 1201 iF iBluck et al.1 1201 11) . iLotz et al.1 JMlT ad- 
dressed the issue of disparate observational merger rates; they 
found that the major reason for these differences is the as- 
sumed timescale in which a merger is observable. Using a 
suite of hydrodynamic merger simulations, they constrained 
the observable timescales of three common merger rate es- 
timators - close galaxy pairs, G/M20, and asymmetry - and 
found that if a physically motivated average observability 
timescale was adopted to calculate the merger rates, then these 
rates become largely consistent. The remaining differences 
between the merger rates are explained by the differences in 
the ranges of mass ratio measured by different techniques and 
the differing parent galaxy selection. 

Additionally, for mass-limited samples (M, > 10'° Mq), 
they found excellent agreement between their merger rates 
from close pairs to several theoretical merger rates. Specif- 
ically, they agreed with the merger rates of iHopkins et al.l 
(I2010al) . who, using a combination of semi-empirical models 
and high-resolution merger simulations, concluded that there 
are enough major mergers, to within a factor of ^ 2, to ac- 
count for the observed growth of the bulge population. They 



argue that previous studies reached different conclusions be- 
cause they assumed incorrect merger timescales, rather, if a 
uniform simulated-calibrated merger timescale is used, then 
many of the se wor ks actually come to their conclusion (see 
also Robaina et al. 2010). 

An alternative way to address whether there are enough 
mergers is to examine the phase in galaxy evolution that is 
predicted to correspond to the period of mergers. Adopting a 
simplified model in which blue star-forming galaxies merge 
and transform into red quiescent galaxies, one would expect a 
short period in which galaxies have intermediate colors, i.e., 
they are in the green valley. A re cent study on the mor pholo- 
gies of green valley galaxies by iMendez et alj (1201 \\) found 
that only 14% of their sample are identified as on-going major 
or minor mergers (using G/M20 and asymmetry parameters), 
which is lower than the 19% merger rate in the blue cloud. 
They further found that most green valley galaxies have disks 
and that 21% have B/T < 0.05, implying that they were un- 
likely to have experienced a recent major merger 

To conclude, while the merger model fits many aspects of 
our data, there are serious open questions, including whether 
major mergers are able to quench, whether there are enough of 
them, and whether they are consistent with the color and mor- 
phologies of green valley galaxies. In the remaining discus- 
sion, we explore other models that may alleviate these prob- 
lems. 

5.2. Disk Instabilities: Violent and Otherwise 

An alternative bulge-building process involves the growth 
of giant c lumps for med via gravit ational instability in gas- 
rich disks (lNoguchiiri 999: Elmegre en et alJl2008tlDekel et alJ 
l2009t ICeverino et al.l 12010 ). These clumps migrate inward 
and coalesce to form a bulge, and simulations suggest that 
galaxies ca n develop class ical bulges with n w 4 during this 
proces^ ( Elmegreen et al. 2008). Recent simulations also 
show that the these clump-origin bulges have central surface 
mass densities comparable to that of our red sequence galax- 
ies (Ceverino et al., in prep). The effectiveness of this in- 
stability, however, is highly dependent on the gas inflow rate 
onto the galaxy (Dekel et al. 2009), which declines with time. 
Thus, we expect this process to be more important at redshifts 
higher (z ^ 2) than that of our sample. 

Although this process may have operated strongly at z 2, 
we stress that this paper concerns a different sample of galax- 
ies at lower reds hift when cond i tions may have changed. 
iBell et al. ( 2004b) and lFaber et 111 (12007) found that the num- 
ber of red sequence galaxies has at least doubled from z ^ 1 to 
z ~ 0, indicating that a fraction of our galaxies at z ~ 0.65 are 
actively migrating to the red seque nce as we view them. Us- 
ing the NEWFIRM survey, Bram mer et all (1201 Ih found that 
the mass density of quiescent galaxies withM* ^3 x 10"^ Mq 
increased by a factor of 10 fro m z 2 to the present day. 
Similarly, IHopkins et alj (1201 Obi) argue that the vast major- 
ity of ellipticals/spheroids do not form through high redshift 
channels. They state that the observed mass density of bulge- 
dominated galaxies at z 2 is only ^ 5% of its z = value, 
and at z ^ 1, is still only ^ 20-35% of its z = value. Thus, 
most bulges are formed at z ^ 1, meaning the majority of our 
red sequence galaxies are recent arrivals. 

Although we expect violent disk instabilities to be increas- 
ingly less frequent at decreasing redshift, owing to lower gas 

However, a recent paper bv llnoue & Saitot] f201 1) claims that these 
clump-oiigin bulges are more akin to pseudobulges in that they exhibit n < 2. 
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fraction, the actual bulge contribution due to this mechanism 
at z < 1 is presently unknown. Hints at intermediate redshift 
suggest that the process may not be limited to high z- For ex- 
ample, Bournaud et al. (201 1) found that half of a sample of 
z ^ 0.7 disk galaxies are clumpy without any merger signa- 
tures, implying that disk instabilities could still be important 
at that redshift. 

Perhaps we need to think more broadly and to recognize 
that the settling of matter to form regular, axisymmetric disks 
is a very lengthy process lasting many billions of years. When 
any non-axisymmetric forces exist, an inevitable outcome is 
that some mass will be driven to the center Furthermore, in a 
general way the higher the degree of non-axisymmetry in the 
potential, the higher the rate of matter inflow will be. At late 
times, non-axisymmetry has become small and the flow rate 
is low, a process that we call "secular evolution", which we 
discuss next. 

Future studies will resolve this question. 

5.3. Secular Evolution 

Secular evolution jKormendv & KennicuttI 120041) consti- 
tutes the weak end of the spectrum of non-axisymmetric pro- 
cesses in disk evolution. This complex of processes involves 
the slow rearrangement of gas (and stellar) mass due to grav- 
itational interactions between the gas and stars within a disk 
galaxy. A variety of relatively weak non-axisymmetric dis- 
turbances, such as bars, ovals, and spiral structure in the 
stellar component, can create non-central gravitational forces 
that add or subtract angular momentum from the gas, which 
responds by moving inward or outward depending on ra- 
dius. The process sweeps inner gas into the center, where 
it forms stars, and pushes gas farther out to larger radii, 
where it can accumulate in tightly wrapped spiral arm s or a 
ring ( Simkin et al.lll980t lKormend y & Kenni cutt 2004). Sep- 
arately, the gas itself may become mildly gravitationally un- 
stable, which raises the velocity dispersion and causes the gas 
to radiate. This net loss of energy must come from some- 
where, and the gas responds by sinking slowl y to the cen- 
ter, in creasing its (negative) potential energy dForbes et alj 
1201 ll) . Overall, these processes push some mass to the cen- 
ter and other mass to the outskirts, thus increasing n. The 
forces are, however, weak and the process is slow, hence 
the term "secular evolution". The total time required would 
be many dynamical t imescales, an d thus at least several Gyr 
jKormendv & Kenni cutt 2004; Fish er et alj|2009l) . 

Although secular evolution may contribute to some of the 
bulge-building taking place at z ^ 0.65, we do not believe 
it is th e major factor. According to iKormendy & Kennicutt 
( 120041) . bars are a major driver of secular evolution at the cur- 
rent epoch. Assuming that bars are also the main drivers at 
higher redshifts, comparing the bar fraction from the past, 
which is 10%-25 % among late-types at z > 0.8 ( Jogee et aP 
I2004t ISheth et all 120081 Herrington et al. in prep), to the 
current epoch, which is 30-60% at z (ISheth et al.ll2008l: 
ICameron et al.l2010t [Masters et al.l20rTI) . implies that secular 
evolution was not a major bulge -building process at z ^ 0.65. 
Additionally, iKoo et al.l (120051) showed that 85% of lumi- 
nous field bulges within this redshift range are red, arguing 
against secular evolution being the dominant bulge-building 
process since they are expected to mainly produce blue bulges 
( iKormendv & Kennicuttll2004l) . 

Finally, the physical structure of bulges built by secular evo- 
lution differs strongly from ones built by mergers, as reviewed 
bv .Kormendv & Kennicutt (.2004.) . The so-called "classical" 



bulges built by mergers resemble small ellipticals embedded 
in disk galaxies. They have high stellar velocity dispersion 
and high vertical extent above the plane, having been "fluffed 
up" by the merger - in other words, they are true spheroids. 
They can also be very massive and contain a considerable 
fraction of the total mass of the galaxy. In contrast, the 
"bulges" built by secular evolution are relatively flat, have ef- 
fective radii of only a few hundred pc , and hav e low fractional 
masse s. Because of these differences. lKormen dv & KennicutS 
('2004') term these structures "pseudobulges". IDrory & Fishen 
(2007) directly compare the properties of these two types of 
bulges. They isolated a sample of nearby, massive, disk galax- 
ies and classified them into classical bulges and pseudobulges 
based on the morphology of the central regions. Confirming 
their disparate nature, they found a clear bimodality in that 
pseudobulge galaxies are much bluer (in the blue cloud or 
green valley), have low central surface brightness, and have 
low global Sersic index (n < 2.5). 

For these reasons, we conclude that the pseudobulges can- 
not be an important contributor to our intermediate redshift, 
high n, red sequence galaxies. They may, however, certainly 
help build the bulges seen in late-type galaxies, becoming in- 
creasingly more important with decreas ing redsh ifts, where 
the bar fraction can be as high as 60% (iShefli et a l. 20080 
In fact, iFi sher & Drorvl ( 12011b show that, by number, 80% 
of the bulges within 1 1 Mpc of the Milky Way are actually 
pseudobulges. But, by mass, they only make up < 10% of 
the total mass de nsity in local spheroids (lAUen et al.l 120061 : 
[Driver et al.ll2007l) . 

In conclusion, we find it helpful to think of the entire fam- 
ily of bulge-building mechanisms as ordered along a "distur- 
bance continuum" from severe to mild, with corresponding 
timescales from short to long and bulge-building rates from 
fast to slow. The members of this continuum consist of ma- 
jor to minor mergers at the strong end, through violent disk 
instabilities, to milder disk instabilities, and finally to weak, 
secular instabilities like bars, spiral arms, and normal star for- 
mation. The dividing line between an externally triggered 
process like mergers and internally triggered dynamical insta- 
bilities is fairly clear, but there is no such dividing line among 
the internal processes - each one shades smoothly into the 
next. However, since every disk is sooner or later subject to 
one or more of these processes, the central densities of disks 
inevitably tend to increase - the only question is, how fast? 

We note that this disturbance continuum is also a contin- 
uum in time, with early galaxies experiencing disturbances 
at the strong end of the spectrum and later galaxies settling 
down to slower, more secular rates. Our AECjIS galaxies ex- 
ist somewhere near the middle of the time continuum, when 
disky galaxies as a class were considerably more disorganized 
and more non-axisymmetric than they are today. This logic 
further supports our conclusion above that secular processes 
were probably not the major bulge-building process in these 
galaxies, with a combination of mergers and stronger disk in- 
stabilities being more likely. However, the exact balance of 
these two processes remains for further study. 

5.4. Morphological Quenching 

The preceding sections focused on mechanisms to increase 
the central densities of disk galaxies, and thus account for one 

This number is c ontroversial, rece nt studies bv lNair & Abrahainl )2O10l) : 
[Masters et alj 1201 lh : lLee et alj 420121) show that the local bar fraction is ^ 
30%. 
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of our major findings, namely, higher central stellar densities 
in quenched galaxies. Even though a variety of stellar build- 
up mechanisms were identified, including mergers and inter- 
nal instabilities, we tacitly assumed regardless of process that 
high density would always lead to the creation of a black hole 
and that feedback from the black hole would quench star for- 
mation. However, the discussion in §5.1 noted a lack of con- 
clusive evidence that AGN feedback actually quenches star- 
formation in disk galaxies. In this section and the next, we 
review two other quenching mechanisms that have been pro- 
posed to operate in central galaxies. 

The first of these is morphological quenching ( Marti g et al.l 
120091) . whereby the steep potential well of a bulge is able to 
drive the Toomre Q parameter above unity and stabilize the 
gas disk against star formation. An attractive aspect of mor- 
phological quenching is that it sets in when n is high, which 
jibes with the structure of quenched galaxies. This mecha- 
nism is unique because it does not require the removal of gas 
or the suppression of the cold-gas supply onto the galaxy. In- 
stead, gas can continue to accrete onto a galaxy yet remain 
inert to star formation owing to the strong central mass con- 
centration. A recent analysis of a set of three high-resolution 
AMR simulations at z ~ 2.3 by Ceverino et al. (2010) demon- 
strates this process. From z — 2.3 to z ^ 1, two of these galax- 
ies are shown to evolve from a gravitationally unstable and 
turbulent disk into a stable system; they attribute this change 
to presence of a dominant stellar bulge. 

An explicit prediction of morphological quenching is that 
red, early-type galaxies could host significant cold gas in 
the amount of a few percent of thek baryonic mass, and 
in fact c omparable to gas fr actions in normal star-forming 
galaxies (" Martig et al.l 120091) . However, this prediction is 
not consistent with observed H 1 properties of nearby early- 
type galaxies. Though H I is frequently detected, espe- 
cially in field galaxies (Morganti et al. 2006; Oosterloo et ^ 
|20l3), the amounts are nearly always low. This is con- 
firmed by the GASS H 1 survey at Arecibo of 1000 
slightly more distant massive galaxies with > 10^" Mq 
dCatinella et alJl20Tot) . which shows that the average H I frac- 
tions of red sequence galaxies are at least ten times lower 
than gal axies of similar mass on the star-forming ma in se- 
quence jSchiminovich et aUlMol; iFabello et alj|20Tlh . The 
same result apphes to molecular H2 in the same galaxies 
jSaintonge et al.l 1201 l l). These recent studies simply con- 
firm what has been known for a long time about early-type 
galaxies, that their absolute gas contents are low comp ared to 
star-formin g galaxies (e.g., vanDriel&van Woerden 119911 
[Roberts & H avnes 1994; Noordermeer et al. 2005). Hence, in 
order to achieve quenching, it is necessary to reduce the frac- 
tional gas content, either by expelling gas or preventing new 
gas from falling in - morphological quenching alone cannot 
do the job. 

For th e nearest galaxies it is possible to map the H I distri- 
bution s ("van Driel & van Woerdenlll9 91l; iNoordermeer et alj 
l2005t Morgant ietaLl 120061; lOosterlo o et all I2010t) . Nor- 
mal lenticular galaxies typically have H I in an outer 
ring with an empty hole in the middle. Several authors 
jvan Driel & van Woerden] 119911; ICortese & Hughes! 120091; 
lOosterloo et al. 2010) have stressed that the surface density of 
gas in thes e rings is well below the critical threshold for sta r 
formation ( lKennicutt|[T989t iSchavdllOOl iBigiel et al.ll2008h . 
which represents the threshold for molecula r H2 formation, 
a tracer of gravitationally unstable gas (e.g. iKrumholz et alj 



I2OI ll) . Moreover, the predicted star formation efficiency at 
such low densitie s is ten times lower t han in normal star- 
forming galaxies (IKrumholz et alj|2012h . which agrees well 
with the low st ar-formation efficiencies seen in the Arecib o 
GASS survey ( Schiminovich et al.ll2010l; IFabello et all201 11) . 
Such low-level star formation has recently been detected in 
the outer disks of normal green-valley SO galaxi es in HST U 
imaging dSalim & Richll2010t ISaUm et al.ll20Tl Fang et al., 
in prep). 

Hence, we reach a very important conclusion about the 
evolutionary track of quenching: low gas content is the un- 
derlying cause of quenching, but star formation shuts down 
even faster as gas content falls below the threshold value, ow- 
ing to the nonlinear relation between cold gas surface den- 
sity and star formation rate. This fall in star-formation effi- 
ciency causes galaxies to redden even faster than expected, 
propelling them rapidly to the red sequenc^3- 

It is fair to point out that all of the data cited above to eval- 
uate morphological quenching comes from nearby galaxies, 
some of which are members of virialized clusters and dense 
groups that are subject to environmental processes such as 
ram-pressure stripping or strangulation. However, many of 
the nearest gal axies are known not to be in clusters (e.g., 
lOosterloo et al.l l2010). and there must also be many field ob- 
jects amongst the thousand or so galaxies in GASS, yet the 
trends are the same. In short, we cannot think of any rea- 
son why morphological quenching would be the key causative 
agent for quenching at higher redshift when it does not appear 
to play that role (even for field galaxies) today. 

5.5. Critical Halo Mass 

The termination of the cold gas supply due to a critical 
halo mass is commo nly referr e d to a s 'halo quenching' 
dSilkl 119771: iR ees & Osti-ikei^ ^Wvf; 'Blumenthal et ajj 
19841; iBirnboim & Dekel 2003 ; Keres et al. 2001 
Dekel & Birnboim 2006; Cattane o et all I2OO6I) . Halos 
below ~ 10'^ Mq are able to accrete gas through cold 
flows while halos above the threshold mass experience a 
virial shock that heats the gas. The hot, diluted gas in 
massive halos is vulnerable to feedback from AGN, which 
effectively halts star formation. An interesting outcome of 
this theory is the ability to generate a hot hal o atmosphere 
that aUows for AGN 'radio m ode' feedback dCroton et al.l 
I2006t [Pekel & Birnboimll2006l) . This mechanism provides a 
way to permanently quench a galaxy, which is desirable since 
the stellar populations of most local ellipticals show no signs 
of rec ent star formation (since at least z ^ 1 ; iThomas et al.l 
120051) . despite continual gas infusion by stars. 

Various studies have explored this theory. For example, 
using a data-driv en, halo-abundance matchin g technique that 
spans < z < 1 . iConrov & Wechsleii (120091) found a gradual 
transition of galaxy properties - including specific star for- 
mation rate - ac ross a halo mass of ~ lO'^ Mq. Recently, 
iMore et al] ( 1201 1 ) inferred the halo mass-stellar mass relation 
based on kinematics of SDSS satellite galaxies and found that 
red, central galaxies, on average, occupy more massive halos 
than blue centrals for fixed luminosity, but show a less ap- 

An interesting corollary comes from tlie fact that UV colors are more 
sensitive to weak star formation than optical colors, and thus galaxies can be 
on the red sequence according toV —B but in the green valley acco rding to 
NUV-r; this is actually seen (e.g.. lSalim & Riclj[2O10l : fSalim eriijlMa) . 
Since galaxies in this paper are classified using U—B, it is possible that some 
of our red sequence objects would appear in the green valley if near-UV color 
were used. 
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preciable difference for fixed stellar mass. iWoo et alJ (l20T2h 
study the quenched fraction of central galaxies as a function of 
halo mass and find that the span of halo masses between 20% 
quenched fraction and 60% quenched fraction is a whole 1.5 
dex. This agrees well with the width of 0.8 dex seen here in 
the overlap region using stellar mass (and the same quenched 
fractions) in Fig. ^ and the theoretical scaling law betwee n 
central and halo mass jKang et alj2005l;ICattaneo et alj2006l) . 
By contrast, the width using our preferred parameter, inner 
surface density, is only 0.5 dex (cf. Fig. 5f). An interpretation 
that emerges from these works is that halo mass has a grad- 
ual and probabilistic effect upon galaxy properties, which is 
therefore very consistent with a wide overlap region in halo 
mass. Instead of a sharp transition at ^ 10^^ Mq, galaxies 
seem to become ready for star-formation quenching around 
this critical value, followed by some type of event that ulti- 
mately triggers quenching. 

This scenario ties in well with our two-stage picture in 
which galaxies "ripen" along the blue cloud, becomi ng m ore 
and more susceptible to quenching as they age. In ^5.11 we 
associated this ripening with decreasing gas content with stel- 
lar mass in the blue cloud, which is also seen in our data as an 
increase in U -B with stellar mass (Fig.|5}5). Lower gas con- 
tent means that less gas needs to be removed in quenching, 
which means in turn that more massive blue cloud galaxies 
are more vulnerable to quenching. Their lower gas content 
stems from the fact that their halo masses are closer to the 
critical value ^ 10'^ Mq, the neighborhood where cold accre- 
tion shuts down. This basic picture does not change if jr^ 
dFranx et alj|200 8) or central velocity dispersion ( Wake et al] 
I2012bh is substituted for stellar mass - any parameter that 
tracks halo mass reasonably well can serve as a ripeness indi- 
cator 

This logic leads to a picture in which the changes in gas 
fraction, color, and star formation rate along the blue cloud 
are caused by the gradual dominance of shock-heated gas over 
cold streams as galaxies near the critical halo mass. A tougher 
challenge, though, is to explain why the actual quenching 
state relates so closely to conditions at the center of the galaxy 
- why is this link so close if the primary governor of cooling 
is out in the halo? 

We have no firm answer to this but offer some speculations. 
Evidently the central conditions either signal, or even trigger, 
a second quenching mechanism and this, plus "natural" halo 
quenching, is what tips a galaxy over the edge. The obvi- 
ous candidate for this second process is AGN feedback, but 
we have stressed that direct evidence for this is still weak. 
It is good to be cautiou s since the ERJS Milky Way simula- 
tion dGuedes et alJ201 lb develops a red bulge and high central 
stellar density quite naturally through mergers and/or internal 
disk evolution. Its star formation rate is falling and it seems 
well on its way to the red sequence, all without benefit of 
AGN feedback. 

Ideal would be some mechanism that both correlates with 
central density and can switch a halo quickly from cold mode 
to hot mode. Some possibilities come to mind. Perhaps AGN 
feedback helps to heat the halo. Perhaps exhaustion of gas at 
the center enhances the ability of stellar winds to sweep gas 
out of the galaxy. Finally, perhaps a merger simultaneously 
builds up central stellar density and triggers a full standing 
shock. Such a transformation is seen in simulations (Dekel 
et al., in prep) where a minor merger triggers a shock that 
expands from the halo center to the virial radius and heats the 



medium. The larger point is that quenching is probably not 
just one factor but a combination of factors that build to some 
critical threshold. 

5.6. The Relationship Between Color and Star Formation 

In this section, we broach the lingering issue of the rela- 
tionship between color and star formation. Throughout this 
paper, we have constantly interchanged these two parame- 
ters, suggesting that color is a good representation of spe- 
cific star formation rate. However, the issue of dust has 
not been thoroughly addressed in our colors. Although we 
have excluded edge-on galaxies (that are presumably highly 
affected by dust) and ensured that our red sequence galax- 
ies are truly quiescent (using the UVJ diagram; see we 
have not actually made any dust correction to our rest-frame 
U-B colors. Therefore, the color trends that we have ex- 
amined throughout this paper may not exactly translate into 
star formation trends. The analysis most affected would be 
our interpretation of how the star formation rates of galax- 
ies behave within the blue cloud. We remarked on the tight 
trend between color and M^/r^ (and M,) within the blue 
cloud in Figs. [Sj) (and |5};) and we proposed that star for- 
mation must decrease with M^,/r^ (or M*). However, this 
color trend may instead be due to more dust in larger galax- 
ies. If this is the case, then a galaxy's star formation would 
be independent of M^jr^ (or M*) within the blue cloud. 
However, work using dust-corrected or dust-robust tnulti- 
wavelength data (e .g., Noeske et al. 2007t JSalim et al.l 120071: 
ISchiniinovich et ak 2007; Zheng et al. 2007) have shown that 
specific star formation rate does indeed decrease with increas- 
ing stellar mass, which is what we have inferred from our 
color-mass diagram (Fig. [S]?). Thus, these works justify our 
subtle assertion that color is a proxy for star formation. Most 
importantly though, our lack of dust-corrected V -B colors 
do not affect our main result that the inner stellar structure of 
galaxies is most related to quiescence, since our quenching 
analysis is based on differentiating galaxies on the blue cloud 
from those on the red sequence, which we have ensured to be 
unaffected by dust (see ®. 

6. CONCLUSION 

In this paper, we analyze a sample of DEEP2/AEGIS galax- 
ies in the redshift range 0.5 < z < 0.8 using H ST /ACS 
V,I images. Our sample has been run through GIM2D, 
a bulge-ndisk decomposition program that gives us in- 
formation on the subcomponents of intermediate red- 
shift galaxies. Using these data, which we provide at 
http: / /people .ucsc. edu/ -echeungl /data .html, our 
goal is to address how quenching depends on galactic struc- 
ture. 

Our methodology is to assess the color correlations of sev- 
eral structural parameters -M*, M*/re, M^/r^, n, and - 
by computing an 'overlap region', which is the band in color- 
parameter space that encompasses both red and blue galax- 
ies. To quantify overlap regions, we calculate the fractional 
number of galaxies within the overlap region AN/N and the 
fractional extent of the region Ax/x; the parameter with the 
smallest overlap region is considered to be the best predictor 
of quiescence. Finding that no parameter is a perfect predic- 
tor of quenching, we explore the overlap region of M*/re for 
secondary color correlations amongst a variety of other pa- 
rameters, including those of bulge and disk. We also consider 
the number of severe outliers from the best predictor of qui- 
escence, n. Our results are the following; 
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1. The Sersic index (n) most sharply discriminates be- 
tween the red sequence and the blue cloud. Eliminating 
dusty, red sequence contaminants ensures that « targets 
quiescent galaxies, not merely red ones. Moreover, the 
color-n diagram resembles a step function, suggesting 
that n is related to a physical quenching threshold. 

2. However, there are exceptions to this general behavior 
- blue galaxies make up « 40% of our n > 2.3 galax- 
ies. Suspecting contamination from starbursts, AGNs, 
and/or other sources of error, we measured central sur- 
face stellar mass densities, which revealed that these 
outliers do not truly belong with the red sequence - 
their Sij^p^ are much lower Central surface density cor- 
rects these outliers and hints that it is the inner structure 
of galaxies that is most related to quiescence. 

3. Red sequence bulges are roughly twice as massive as 
blue cloud bulges at the same galaxy stellar mass (and 
M^,/re), yet also roughly twice as small. This structural 
difference shows that most blue galaxies at the observed 
epoch do not simply fade onto the red sequence. Rather, 
the high values of n and Sij^pj, on the red sequence must 
reflect a net migration of existing stars toward the center 
of the system or the formation of new stars at the center 
This restructuring either causes quenching itself or is 
closely related to the process that does. 

4. While in the blue cloud before quenching, a galaxy's 
star formation rate is most closely correlated with 
M*/re (or with M*). 

These results suggest that galaxies evolve toward the red se- 
quence in a two-stage process. In stage one, a galaxy is star- 
forming in the blue cloud at a rate that correlates with global 
parameters like M^/rf^ (or M*). Since these in turn reflect 
halo mass, a broad conclusion might be that star-formation is 
controlled by the galaxy's halo while stars are still forming. 
As stellar mass increases, the halo mass also increases, ulti- 
mately approaching the critical value ^ 10^^ M©, where cold 
flows begin to give way to hot gas, which cannot accrete. The 
gas content of the galaxy begins to fall as it nears this critical 
value, and colors redden. 

The result of stage one is a galaxy that is increasingly vul- 
nerable as time goes by to the onset of a second quenching 
process. This second process must be closely associated with 
bulge-building, and central stellar density, Sji^p;,, must in- 
crease. AGN feedback is an obvious candidate for this second 
process, but direct evidence for this is still weak. Also unclear 



is the exact mechanism that leads to the central mass build-up 
at the center of the galaxy. We have stressed that galaxies, par- 
ticularly at high redshift, are far from axially symmetric and 
that any non-axisymmetry leads inevitably to an exchange of 
angular momentum and/or loss of energy, which causes some 
stars and gas to move inward. Major mergers sit at the strong 
end of this "disturbance continuum", secular evolution pro- 
cesses sit at the other end, and minor mergers, violent disk 
instabilities, and milder disk instabilities sit in the middle. 
Given that our galaxies lie at z ^ 0.65, where galaxies are 
still moderately disturbed, it is unlikely that secular evolution 
plays a major role in them. More likely is some combination 
of mergers and disk instabilities, which collectively are proba- 
bly strong enough and frequent enough to do the job. Mergers 
in particular have several well known advantages: they build 
bulges at the right stellar mass, they naturally build a con- 
centrated stellar spheroid, they drive a lot of gas quickly to 
the center that can power an AGN or starburst, and they may 
be able to quickly switch a marginally cooling halo into hot 
mode. A problem for major mergers is the large fraction of 
galaxies in the green valley that are disky, but these might be 
explained by the other mechanisms or they might be reddened 
members of the blue cloud. 

One conclusion seems clear, and that is that moving into the 
green valley and thence to the red sequence requires a lower- 
ing of fractional cold gas content. This can only be achieved 
either by expelling gas or by preventing its infall. Exactly 
how this happens is still not clear, but at least some of the 
parameters surrounding the process are better known. 
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APPENDIX 
SURFACE BRIGHTNESS LIMITS 

One might be concerned that our 'final' sample may be missing low surface brightness galaxies; in this Appendix, we show 
that this is not the case. We demonstrate our 'final' and 'starting' samples' surface brightness limits in Fig.[TO]by plotting the 
effective radius in arcseconds vs. V and / (see Ofor sample definitions). The master GIM2D sample is plotted in gray and the 
'starting' and 'final' sample are overplotted in red and blue, respectively. In each panel, we draw two lines of constant surface 
brightness to give an idea of our samples' surface brightness boundaries. There is obviously an edge to the distribution of the 
master GIM2D sample at low surface brightnesses. Whether this is a selection limit or a reflection of where real galaxies lie is 
unclear But it is clear that the master GIM2D sample contains galaxies with the lowest surface brightnesses, indicating that our 
'starting' and 'final' samples are sensitive to them. 

In these plots, the relevant difference between the 'starting' and 'final' sample is the imposed mass limit of the 'final' sample. 
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Figure 10. Semimajor axis effective radius in arcsecond is plotted against V and /. The two lines of constant surface briglitness in eacli panel roughly marks the 
edge of the 'starting' and 'final' samples' surface brightness distribution. The contours represent the number density of the inaster GIM2D sainple. 



The main concern in terms of sample completeness boils down to whether the 'final' sample is missing massive low surface 
brightness galaxies. It is clear that there are many galaxies in 'starting' sample (and the master G1M2D sample) with much lower 
surface brightnesses than those in the 'final' sample. Hence if such massive low surface brightness galaxies existed, the 'final' 
sample should certainly contain them, thus our 'final' sample does not miss low surface brightness galaxies. 

GIM2D MEASUREMENT QUAILITY 

In this appendix, we show that the errors on the GIM2D measurements are well-behaved even at the sample limits. The most 
important quantities are the GIM2D model HST/ ACS V and / magnitudes since every structural measurement is based on them. 
Fig. [TTIplots the V and / 1-<t errors as a function of V, I, and log M*. These errors are confidence limits derived through full 
Monte Carlo propagations of the parameter probability distributions computed by GIM2D (Simard et al. 2002). 

Shown are the master GIM2D catalog (gray), 'starting' sample (red), 'final' sample (blue), and galaxies within the 'final' 
sample with Sersic indices n < 1 (pink). The 5-(t limiting magnitudes for the HST/ACS images are V = 26.23 and / = 25.61, 
consistent with the master GIM2D distribution. Within our 'final' sample, the limits of the magnitude distributions are w 0.5 
mag and « 1 .0 mag brighter than those limits, respectively. And most importantly, the V and / magnitude errors at the mass and 
magnitude limit of the 'final' sample are small, < 0.06 mag, with the majority < 0.03 mag. 

To further illustrate the quality of the 'final' sample's GIM2D fits, we plot in Fig.[T2]the 1-cr error of the Sersic index, fractional 
1-a error of the effective radius, and the l-a error of the bulge-to-total ratio B/T (in the / band) as a function of V, I, and M*. 
The errors of the Sersic index measurements are almost all < 0.2. A typical late-type galaxy has n ^ I while a typical early- 
type galaxy has n ^ A. An error of 0.2 on n will not affect this division and hence these errors are tolerable. The errors on the 
effective radius measurements are presented as fractional errors, i.e., l-cr effective radius error divided by the effective radius. 
Note that these radii are semimajor axis effective radii and not circularized effective radii. The fractional errors are small, almost 
all are < 6%. The errors on the B/T are almost all < 0.05. Most star-forming, blue galaxies in our 'starting' sample, which 
are presumably late-type, have B/T < 0.10 while most quiescent, red galaxies in our 'starting' sample, which are presumably 
early-type, have B/T > 0.40. Thus an error of 0.05 will not push blue galaxies to B/T values of red galaxies and vice versa. In 
summary, these plots demonstrate that our 'final' sample contain quality GIM2D measurements down to the mass and magnitude 
limits of the 'final' sample. 

SERSIC INDEX BIAS? 

One might be concerned that the Sersic index may be sensitive to surface brightness and stellar mass. Fig.[T3]plots the Sersic 
index against these quantities. The top row is the fi inV and /. The vertical green dot-dashed line represents the approximate 
edge of the 'final' sample's surface brightness distribution as seen in Fig. [TO] The contours represent the number density of the 
master GIM2D catalog, which clearly shows that the majority of these galaxies have low n. At faint /i, the 'starting' sample 
clearly has less n > 3 galaxies compared to the n < 3 population, but we would argue that is a reflection of the universe, i.e., high 
n galaxies generally have bright fi. For the 'final' sample, however, this proportional dearth of n > 3, faint /i galaxies seen in the 
'starting' sample seems to have been greatly reduced. 

We also plot the Sersic index against stellar mass in the lower panel. At low stellar masses, the 'starting' sample also has 
less high n galaxies than low n galaxies. However, this absence does not seem as severe as with faint fi. And within our 'final' 
sample, it seems that this disparity in Sersic populations have decreased. Moreover, it is important to note that even at the lowest 
masses, our 'final' sample still contains Sersic index values across the entire possible range, indicating that GIM2D does not 
preclude high n fits for low mass galaxies. In fact, the absolute count of galaxies with Sersic index n > 3 in every mass range is 
approximately the same. Hence there is no predisposition against n > 3 galaxies in our 'final' sample. The fact that our 'final' 
sample has a large concentration of low n galaxies is most Ukely due to the fact that there are many late-type galaxies at low 
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masses. 

One might wonder how the low n galaxies fared near the limits of the data. As can be seen in Fig.[T3] the low n galaxies make 
up a significant fraction of our 'final' sample. To see the reliability of the measurements of these galaxies, we've highlighted this 
population (« < 1 galaxies in the 'final' sample) in magenta in several of the previous figures, specifically. Fig. [TT] and [12] In 
most cases, these low n galaxies have errors consistent with the rest of the population, i.e., low n galaxies fare fine near the limits 
of our 'final' sample. 



REFERENCES 



Abraliam, R. G., Valdes, R, Yee, H. K. C, & van den Bergli, S. 1994, ApJ, 432, 75 

Alexander, D. M., Swinbank, A. M., Smail, I., McDermid, R., & Nesvadba, N. P. H. 2010, MNRAS, 402, 2211 

Allen, P. D., Driver, S. P, Graham, A. W., et al. 2006, MNRAS, 371, 2 

Ammons, S. M., Rosario, D. J. V., Koo, D. C, et al. 2011, ApJ, 740, 3 

Amouts, S., Walcher, C. J., Le Fevre, O., et al. 2007, A&A, 476, 137 

Baldwin, J. A., Phillips, M. M. & Terlevich, R. 1981, PASP 93, 5 

Barmby, P, Huang, J.-S., Ashby, M. L. N., et al. 2008, ApJS, 177, 431 

Barro, G., Perez-Gonzalez, R G., Gallego, J., et al. 2011a, ApJS, 193, 13 

Barro, G., Perez-Gonzalez, R G., Gallego, J., et al. 2011b, ApJS, 193, 30 

Bell, E. F, & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. F, Mcintosh, D. H., Barden, M., et al. 2004a, ApJ, 600, LI 1 

Bell, E. F, Wolf, C, Meisenheimer, K., et al. 2004b, ApJ, 608, 752 

Bell, E. F, Phleps, S., Somerville, R. S., et al. 2006, ApJ, 652, 270 

Bell, E. F 2008, ApJ, 682, 355 

Bell, E. F, van der Wei, A., Papovich, C, et al. 2012, ApJ, 753, 167 
Bertin, E., & Amouts, S. 1996, A&AS, 117, 393 

Bezanson, R., van Dokkum, P G., Franx, M., et al. 2011, ApJ, 737, L31 
Bigiel, F., Leroy, A., Walter, R, et al. 2008, AJ, 136, 2846 
Bimboim, Y., & Dekel, A. 2003, MNRAS, 345, 349 



The Dependence of Quenching upon Inner Galactic Structure 




16 19 20 21 22 25 24 25 26 27 2B 29 
ACS V mognrtude 




18 19 20 31 32 23 34 35 36 37 36 39 
ACS V mognrtude 




IS 19 2 21 22 2 3 24 2 5 26 2 7 2 8 29 
ACS I magnitude 



0.6 r 

0.5 \ 



? 0.4 r 




log M. 

0,1 Of ' • — ' — 



0.08 ■ 




log M. 




Bimboim, Y., & Dekel, A. 2011, MNRAS, 415, 2566 

Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 594, 186 

Blanton, M. R., & Rowels, S. 2007, AJ, 133, 734 

Bluck, A. F. L., Conselice, C. J., Buitrago, R, et al. 2011, arXiv:l 1 1 1.5662 
Blumenthal, G. R., Faber, S. M., Primack, J. R., & Rees, M. J. 1984, Nature, 311, 517 
Borch, A., et al. 2006, A&A, 453, 869 

Boumaud, F, Jog, C. J., & Combes, F 2007, A&A, 476, 1 179 
Boumaud, F, Juneau, S., Le Floc'h, E., et al. 201 1, arXiv: 1 1 1 1.0987 
Brammer, G. B., Whitaker, K. E., van Dokkum, R G., et al. 2009, ApJ, 706, LI 73 
Brammer, G. B., Whitaker, K. E., van Dokkum, R G., et al. 201 1, ApJ, 739, 24 
Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 
Brown, M. J. 1., Dey, A., Jannuzi, B. T., et al. 2007, ApJ, 654, 858 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bundy, K.,et al. 2006, ApJ, 651, 120 

Bundy, K., Fukugita, M., Ellis, R. S., et al. 2009, ApJ, 697, 1369 
Bundy, K.,et al. 2010, ApJ, 719, 1969 

Cano-Diaz, M., Maiolino, R., Marconi, A., et al. 2011, arXiv:1112.3071 

Cameron, E., Carollo, C. M., Oesch, R, et al. 2010, MNRAS, 409, 346 

Capaccioli, M. 1989, World of Galaxies (Le Monde des Galaxies), 208 

Catinella, B., Schiminovich, D., Kauffmann, G., et al. 2010, MNRAS, 403, 683 

Cattaneo, A., Dekel, A., Devriendt, J., Guiderdoni, B., & Blaizot, J. 2006, MNRAS, 370, 1651 

Ceverino, D., Dekel, A., & Boumaud, F 2010, MNRAS, 404, 2151 

Ciotti, L., Ostriker, J. R, & Proga, D. 2009, ApJ, 699, 89 

Coil, A. L., Newman, J. A., Kaiser, N., Davis, M., Ma, C.-R, Kocevski, D. D., & Koo, D. C. 2004, ApJ, 617, 765 
Coil, A. L., Weiner, B. J., Holz, D. E., Cooper, M. C, Yan, R., & Aird, J. 2011, arXivil 104.0681 
Collister, A. A., & Lahav, O. 2004, PASP, 116, 345 
Conroy, C, & Wechsler, R. H. 2009, ApJ, 696, 620 
Conselice, C. J. 2006, ApJ, 638, 686 

Conselice, C. J., Bundy, K., U, V., et al. 2008, MNRAS, 383, 1366 
Cooper, M. C, Aird, J. A., Coil, A. L., et al. 201 1, ApJS, 193, 14 



24 



Cheung et al. 



Final Somple 
Starting Sampie 

MASTER GIM2D 

6N% 



■D 

^ 4 

u 

3 




Finai Sompie 
Storting Sample 

MASTER GIM2D 



V 

c 
u 

I- 

u 




16 IB 20 22 24 26 
/I (V) [mag arcsecond"^] 



28 




16 18 20 22 24 26 
fi (l) [mag arcsecond"^] 



28 



Figure 13. Sersic index is plotted against ^i(V) and in the top row. Contours represent the number density of the GIM2D m aste r sample in gray. The vertical 
green dot-dashed line represents the approximate edge of our 'final' sample's surface brightnesses distribution (as seen in Fig. llOt . The bottom panel plots the 
Sersic index against stellar mass. The 'final' sample does not show a strong disparity of high and low n populations at faint (low) surface brightness (mass). 
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Table 1 

Galaxy Properties 





Q 


z 


Mb 


Mb 


U —n 


U — D 












/c-correct 


t,qns.|i]|4] 


^-correct 


t,qn.|2J 


J. Huang et al., in prep 


TJnnc- 1/^11111/1 1 
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Note. — Twenty randomly-selected galaxies from the catalog is shown. The entire catalog is available at: 
]http : //people . UCSC.edu/~echeungl/ data, html Col. (l): Unique DEEPID. Col. (2): Spectroscopic 
redshift quality code. Values of 3 and above are secure spectroscopic redshifts {see Newman et al. 2012 for more details) and they 
comprise the spectroscopic sample. For values below 3, we use the photometric sample and their coixesponding photometric redshifts. 
Col. (3): Best available redshift For those with (2 > 3, these are spectroscopic z, those with Q <2 are photometric z- Col. (4): 
Absolute fi-band magnitude derived from ^:-correct. About 7% have large eiTor measurements, we use vales from Col. (5) for these 
objects. Col. (5): Absolute S-band magnitude derived from Eqns.[3]|4] Col. (6): U -B rest-frame color derived from A:-coiTect. About 
^ 7% have large enor measurements, we use vales from Col. (7) for these objects. Col. (7): U — B rest-frame color derived from 
Eqn.[5] Col. (8): Stellar mass from J. Huang et al., in prep. 10% of our sample have no conesponding \ they are marked by 
-1.00. We use Col. (9) for these objects. Col. (9): Stellar mass derived from Eqns.[6][3]|4] 



Table 2 

Subcomponent Properties 
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Note. — The same twenty randomly-selected galaxies 

from Table fll is shown. T he entire ca talog is available at: 

[http : /7peopie .ucsc.edu/~echeungl/ data, htmll 

Col. (1): Unique DEEPID. Col. (2): Absolute B-band magnitude of bulge 
derived from Egn.llBI Col. (3): Absolute B-band magnitude of disk derived 
from Eqn. |3|4| Col. (4): U —B rest-frame color of bulge derived from Eqn.fS] 
Col. (5): U — B rest-frame color of disk derived from Eqn.fs] Col. (6): of 
bulge derived from Eqn.[3][6] Col. (7): of disk derived from Eqn. [3][6| 
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Table 3 

GIM2D: Single n Catalog 
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Note. — The same twenty randomly-selected galaxies 

from Table [Tl is shown. The entire catalog is available at: 

|http : / /people . ucsc. edu/~echeungl / data . html Eveiy 
value in this table is from the GIM2D decomposition that only fits for a single Sersic 
index. These are the primary data we use throughout this paper, including the basic V 
and / of the galaxy. Col. (1): DEEPID. Col. (2): \/-band magnitude of galaxy. Col. (3): 
/-band magnitude of galaxy. Col. (4): Sersic index of galaxy. Col. (5): EUipticity of 
galaxy, e = \-bja. b = semiminoraxis. a = semimajoraxis) Col. (6): Effective 
radius of the galaxy measured along the major axis in units of pixels. Col. (7): x~ of fit 
the V-band. Col. (8): of fit in the /-band. 



Table 4 

GIM2D: n = 4 Bulge Catalog 
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23.07 


22.47 


22.71 


0.46 


0.56 


19.41 


5.88 


0.55 


0.58 


1.13 


1.05 


12015606 


24.46 


22.67 


22.88 


21.95 


0.16 


0.30 


11.22 


10.18 


0.48 


0.44 


1.08 


0.97 


13058131 


24.43 


23.65 


22.59 


22.42 


0.33 


0.46 


7.95 


19.11 


0.57 


0.55 


0.97 


0.88 


12016156 


28.30 


23.97 


25.40 


23.57 


0.02 


0.16 


18.57 


15.65 


0.40 


0.30 


1.11 


1.00 


13012297 


23.47 


25.07 


21.80 


23.52 


0.81 


0.83 


16.04 


6.93 


0.55 


0.58 


0.89 


0.82 


12007918 


23.45 


23.19 


22.35 


22.74 


0.44 


0.59 


19.23 


5.87 


0.55 


0.56 


1.16 


1.07 


13011795 


24.78 


25.27 


23.12 


24.09 


0.61 


0.71 


6.49 


12.72 


0.66 


0.65 


0.87 


0.81 


13064645 


24.95 


25.15 


23.56 


23.36 


0.55 


0.45 


2.26 


10.60 


0.69 


0.71 


0.94 


0.88 


12016799 


30.56 


23.11 


28.18 


21.85 


0.00 


0.00 


13.39 


11.62 


0.36 


0.38 


1.17 


1.14 


12023870 


23.39 


24.42 


22.83 


24.18 


0.72 


0.78 


21.91 


1.89 


0.66 


0.71 


1.14 


1.15 


12004470 


23.32 


23.27 


22.90 


22.94 


0.49 


0.51 


11.73 


3.03 


0.73 


0.74 


1.53 


1.42 


12020439 


24.77 


24.01 


23.65 


23.61 


0.33 


0.49 


13.89 


3.54 


0.65 


0.72 


1.02 


1.06 



Note. — The same twenty randomly- selected galaxies _ from Table is shown. The entire catalog is available at: 
fhttp :/ /people . UCSC . edu/~echeungl /data . html,| Every value in this table is from the GIM2D decomposition 
using a bulge and disk with Sersic index of « = 4 and n = \, respectively. Col. (1): DEEPID. Col. (2) & (3): V-band magnitude of bulge and disk. 
Col. (4) & (5): /-band magnitude of bulge and disk. Col. (6): B/T in V-band. Col. (7): B/T in /-band. Col. (8): Effective radius of major axis of 
bulge in units of pixels. Col. (9): Scale length of disk measured along major axis in units of pixel. Col. (10): Concentration with a. = 0.3 in V. Col. 
(11): Concentration with a = 0.3 in /. Col. (12): Chi-squared of fit inV. Col. (13): Chi-squared of fit in/. 
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Table 5 

GIM2D: n = 2 Bulge Catalog 



DEEPID 


V 


/ 


V 


/ 


B/T 


B/T 






C 


C 








Bulge 


Bulge 


Disk 


Disk 


Bulge 


Disk 
















V-band 


/-band 


(pixels) 


(pixels) 


y-band 


/-band 


y-band 


/-band 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


13049654 


23.33 


22.38 


22.51 


21.90 


0.29 


0.36 


24.61 


10.97 


0.47 


0.47 


1.11 


1.00 


13018599 


25.53 


24.55 


23.41 


23.23 


0.29 


0.46 


20.65 


8.46 


0.48 


0.58 


1.05 


1.12 


12020067 


23.02 


23.16 


22.43 


22.70 


0.53 


0.56 


23.38 


4.43 


0.53 


0.55 


1.05 


0.98 


12007757 


24.66 


23.79 


23.48 


23.18 


0.31 


0.43 


14.76 


35.53 


0.35 


0.36 


1.04 


0.98 


13040619 


24.40 


23.95 


23.20 


22.18 


0.40 


0.28 


1.92 


9.09 


0.62 


0.69 


1.01 


0.90 


13048556 


22.99 


25.33 


22.59 


23.80 


0.90 


0.75 


28.14 


30.87 


0.50 


0.48 


1.11 


0.96 


13049852 


23.48 


22.48 


22.14 


21.05 


0.28 


0.27 


3.13 


15.79 


0.61 


0.64 


0.93 


0.88 


13026131 


24.93 


24.34 


23.73 


23.68 


0.37 


0.49 


19.07 


9.98 


0.41 


0.46 


0.95 


0.93 


12020067 


23.02 


23.16 


22.43 


22.70 


0.53 


0.56 


23.38 


4.43 


0.53 


0.55 


1.05 


0.98 


12015606 


23.71 


22.89 


22.36 


22.30 


0.32 


0.49 


10.64 


11.16 


0.48 


0.44 


1.07 


0.97 


13058131 


24.81 


23.55 


22.95 


22.23 


0.24 


0.34 


4.94 


17.80 


0.57 


0.55 


0.97 


0.88 


12016156 


24.67 


24.31 


23.93 


23.91 


0.42 


0.49 


27.86 


30.96 


0.40 


0.30 


1.11 


0.96 


13012297 


23.44 


-99.99 


21.80 


-99.99 


1.00 


1.00 


10.84 


6.88 


0.55 


0.58 


0.91 


0.88 


12007918 


23.33 


23.54 


22.77 


22.68 


0.55 


0.48 


9.21 


7.01 


0.55 


0.56 


1.17 


1.08 


13011795 


25.30 


24.90 


23.63 


23.50 


0.41 


0.47 


3.11 


13.02 


0.66 


0.65 


0.87 


0.80 


13064645 


25.29 


24.90 


23.93 


23.21 


0.41 


0.34 


1.36 


8.69 


0.69 


0.71 


0.94 


0.89 


12016799 


23.41 


23.85 


21.90 


22.78 


0.60 


0.69 


44.45 


9.40 


0.36 


0.38 


1.11 


1.06 


12023870 


23.71 


24.13 


23.12 


23.82 


0.60 


0.66 


21.50 


1.75 


0.66 


0.71 


1.14 


1.15 


12004470 


23.50 


23.19 


23.00 


22.89 


0.43 


0.47 


14.01 


2.61 


0.73 


0.74 


1.48 


1.37 


12020439 


24.56 


24.16 


23.63 


23.72 


0.41 


0.52 


12.21 


3.11 


0.65 


0.72 


0.98 


1.05 



Note. — The sam e twenty randomly- selected galaxies from Table \T] is shown. The entire catalog is available at: 

|http : / /people . UCSC.edu/-echeungl / data, html. Eveiy value in this table is from the GIM2D decomposition 
using a bulge and^isk with Sersic index of n = 2 and n= I, respectively. Col. (1): DEEPID. Col. (2) & (3): V-band magnitude of bulge and disk. Col. 
(4) & (5): /-band magnitude of bulge and disk. Col. (6): B/T in V-band. Col. (7): B/T in /-band. Col. (8): Effective radius of major axis of bulge 
in units of pixels. Col. (9): Scale length of disk measured along major axis in units of pixel. Col. (10): Concentration with a = 0.3 in V. Col. (11): 
Concentration with a = 0.3 in /. Col. (12); Chi-squared of fit in V. Col. (13): Chi-squared of fit in /. 



